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INTEGRAL CONVEX POLYHEDRA 

AND 

AN APPROACH TO INTEGRALIZATION* 

Abstract 



Many combinatorial optimization problems may be formulated as 
integer linear programming problems — that is, problems of the form: 
given a convex polyhedron P contained in the non-negative orthant 
of n-dimensional space, find an integer point in P which maximizes 
(or minimizes) a given linear objective function. Well known 
linear programming methods would suffice to solve such a problem if: 

(i) P is an integral convex polyhedron, or 

(ii) P is transformed into the integral convex polyhedron 
that is the convex hull of the set of integer points 
in P, a process which is called integralization. 

This thesis provides some theoretical results concerning integral 
convex polyhedra and the process of integralization. Necessary and 
sufficient conditions for a convex polyhedron P to have the integral 
property are derived in terms of the system of linear inequalities 
defining P. A number- theoretic method for integralizing two-dimensional 
convex polyhedra is developed which makes use of a generalization of 
the division theorem for integers. The method is applicable to a 
restricted class of higher dimensional polyhedra as well. 



*This report reproduces a thesis of the same title submitted to the 
Electrical Engineering Department, Massachusetts Institute of Technology, 
in partial fulfillment of the requirements for the degree of 
Doctor of Philosophy, August 1970. 
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CHAPTER 1 



INTRODUCTION 



It is well known that many optimization problems of a combinatorial 
nature may be formulated as integer programming problems. Several authors 
have written on this subject, including Balinski [1] and [2], Dantzig [5], 
Hoffman [14], and Hu [17]. Outstanding examples are the problems of 
constructing a maximal flow in a transport network and constructing a 
maximal matching in a bipartite graph. In fact, these two examples lead 
very naturally to integer linear programming formulations which may be 
solved by standard linear programming methods. Other combinatorial 
optimization problems, such as the traveling salesman problem or the 
problem of properly coloring the regions of a map with a minimum number 
of different colors, lead to integer linear programming formulations 
which cannot be solved by standard linear programming methods. They 
require different and more complex solution methods. 

In 1958, Gomory [8] developed the first general algorithm for 
solving integer linear programming problems. There has been considerable 
interest in integer linear programming since Gomory 's pioneering work, 
and a variety of methods for dealing with the problem have been 
developed. Balinski [1] and [2] presents an extensive survey of many 
of these methods and describes some computational experience with them. 
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In 1965, Gomory [9] brought a new algebraic approach to bear on 
the problem. He showed that when certain simplifying assumptions are 
made, the integer linear programming problem can be transformed into 
an optimization problem involving a finite Abelian group. The 
algebraic ideas contained in his approach have been further developed 
and applied by Gomory [10] and [11], Shapiro [23] and [24], Hu [16], 
and others. 

Although a good deal of success in dealing with integer linear 
programming problems has been achieved, we believe there remains 
much to be learned in this area. Motivated primarily by the 
importance of integer linear programming as a tool for solving 
difficult combinatorial optimization problems, this thesis seeks 
to provide some relevant theoretical results which may contribute 
to a better understanding of the problem. 

The thesis is organized into seven chapters. In the remainder 
of this first chapter we review some n-dimensional geometric concepts 
and definitions related to convex polyhedra. We then give some 
graphical examples to illustrate the concept of an integral convex 
polyhedron and the concept of transforming a convex polyhedron into 
an integral convex polyhedron — a process which we call integralization. 
We conclude the chapter with a discussion of the relevance of integral 
convex polyhedra and the process of integralization to the integer 
linear programming problem, and a description of some of the related 
research of others. 



In Chapter 2 we derive several results concerning convex polyhedra 
having the integral property. In particular, we present a necessary 
and sufficient condition for a system of linear equations with integer 
coefficients and integer right hand sides to have an integer solution. 
We then use this to extend the results of Hoffman and Kruskal [15] 
and obtain a necessary and sufficient condition for a particular 
convex polyhedron to have the integral property. 

In Chapter 3 we define a new kind of integral property which we 
call the totally integral property. This definition is motivated by 
two considerations. First, the class of totally integral polyhedra 
is a proper subclass of the class of integral polyhedra. Second, 
the class of linear programming polyhedra having the totally integral 
property can be characterized algebraically, by applying some of the 
group theoretic ideas developed by Gomory [9]. 

In the next three chapters we turn to the process of integralization. 
We shall classify polyhedra using two parameters n and k. By an 
(n,k)-polyhedron we shall mean an n-dimensional polyhedron defined 
by a system of linear inequalities whose coefficient matrix has rank k, 
where k ^ n. 

In Chapter 4 we show that an (n,k)-polyhedron, where n is 
arbitrary, n s k, can always be transformed into a (k,k) -polyhedron 
by a transformation which preserves integralization. This allows us to 
focus on the problem of integralizing (k,k) -polyhedra without loss 
of generality. As we shall see, if we know how to integralize a 
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(k,k) -polyhedron then we know how to integralize an (n,k) -polyhedron 
for arbitrary n ;> k. We then discuss the simple problem of 
transforming and integralizing (n,l)-polyhedra. We show that this 
is accomplished by using the greatest integer function applied to 
rational numbers, an application of the division theorem for integers. 

In Chapters 5 and 6 we deal with integralization of (2,2)-polyhedra. 
We introduce a generalization of the division theorem for integers 
which applies to ordered pairs of integers and use this, together 
with an associated division process, to integralize such polyhedra. 

Chapter 7 concludes the thesis with a discussion of the difficulty 
we encountered in attempting to generalize the integralization process 
to make it apply to (k,k) -polyhedra, for k >3. We then offer some 
suggestions for future research in this area. 



1.1 Geometric Concepts and Definitions 

In this section we review some n-dimensional geometric concepts 
and definitions. Most of the material is well known, with the exception 
of the subsection on integral convex sets, and may be found in one or 
more of the following references: Benson [3], Bonnice and Klee [4], 
Goldman [6], Goldman and Tucker [7], Grunbaum [12], Hadley [13], 
Hoffman and Kruskal [15], Klee [19], and Weyl [26]. 



Let R be the real field. Let R be the n-dimensional vector space 
over R. We refer to x £ R , x = (x.. , x„, . . . , x ) , as a real point 
or real ii-vector (column). 



1.1.1 Dimension. Boundedness. and Convexity 

01 k U 1 ,. 1 *. 

Let x , x , ..., x be real n-vectors. A real n-vector x is 

1 k 
said to be a linear combination of x , x , . . . , x if 

X = X n X + A., X + « • • + A. X 

1 k 

for some real numbers A.„, A. . .., \. . If the A.'s satisfy the 
constraint, 



1 k 
then x is said to be an affine combination of x , x , ..., x . 

If the A.'s satisfy the additional constraint, 



A.. :> 0, for i = 0, l,...,k 

. ^ , ^..^01 k 

then x is said to be a convex combination of x , x , . . . , x . 

1 k 
A set X = {x , x , . . . , x } is said to be linearly dependent 

1k 
if 0, the all zero n-vector, is a linear combination of x , x , . . . , x 

in which some A.. ^ 0. If, in addition, 
l ' 
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then X is said to be af finely dependent . If X fails to be linearly 
(affinely) dependent then X is said to be linearly ( affinely) 
independent . We shall say that X is a dependent ( independent ) 
set iff X is affinely dependent (independent). It can be 

shown that X is a dependent (independent) set iff the set 

r 1 ° 2 ° k 0, . ,. , J J ,. , , 

[x -x , x -x , . . . , x -x ] is linearly dependent (independent). 

3 

Example 1-1 Independent sets in R are: a single real point, 

two distinct real points, three real points not contained in a 
single line, four real points not contained in a single plane. 



□ 



For the following definitions we let X be any subset of R , 
either a finite subset or an infinite subset. 

X is k-dimensional . where <, k <, n, if X contains an independent 
set of k+1 real points but none of k+2 real points. 

X is bounded if there exist two real n-vectors x and x such 
that, for all x £ X, x ^ x ^ x , where <. is componentwise. Otherwise 

X is unbounded . 

1 2 
X is convex if, for all x , x £ X and for all real \, <, X ^ 1, 

X x + (1-X) x £ X. In words, a convex subset X s R is one which 

contains the line segment joining any two of its points. 

An intersection of convex subsets of R is again a convex subset 

of R . 
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■J 

Example 1-2 In R , any single real point is a O-dimensional bounded 

convex set. R itself is a 2-dimensional unbounded convex set. Fig. 1-1 

2 
shows several subsets of R which have various combinations of these 

properties. |— . 



Bounded 



Unbounded 



-v /- 






> 1-dimensional 



> 2-dimensional 



V 



— — y— — 

Convex 



Non-convex 



J 



ii.%. i-i 



12 



1.1.2 Convex Hull 

Given an arbitrary subset X c R ? which may or may not be convex, 
it is always possible to find a convex subset of R which contains X. 
R itself will do the job. The convex hull or convex closure of X, 
denoted by H(X), is defined to be the intersection of all convex 
subsets of R which contains X. H(X) is unique. Intuitively, H(X) 
is the smallest convex subset of R which contains X. If X is 
itself convex, then H(X) = X. 



Example 1-3 Fig. 1~2 illustrates several subsets of R and, 
immediately to the right of each, their convex hulls. 



□ 



(a) 





(d) 



(b) 






(e) 



(c) 



Fig. 1-2 
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If X is a finite subset of R n , X = {x , x , . . . , x™} , then it can 
be shown (see Benson [3]) that the convex hull of X is given by, 



m m 



H(X) = { x g R n f x = 7 X x 1 , \ > 0, Y \ = 1 } (1-1) 

i=0 i=l 

where \q, \^, . .., X are real numbers. In words, (1-1) says that 
H(X) is the set of all convex combinations of X. 

Another characterization of the convex hull of a subset of R n , 
similar to (1-1) but more general in that it applies to infinite 
subsets as well, is given by Caratheodory's Theorem (see Bonnice 
and Klee [4]). Given a subset X c R n , the theorem states that 



n n 

H(X) = ■[ x € R n ( x = T \ ± x 1 , \. :> 0, y \ ± = 1, x 1 € x} (1-2) 

i=0 i=0 

where X Q , \ ..., \ are real numbers. In words, (1-2) states that 
H(X) is the set of all convex combinations of all finite subsets of 
X containing n+1 or fewer points. 



1.1.3 Flats 

Let X = {x , x , . . . , x } be an independent set of real points 
in R , where <. k <; n. A k-flat n in R is a k-dimensional convex 
subset of R , 
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T = { X 6 ^ I K = J \ ± x\ 7 \ = 1 } 



(1-3) 
1=0 i=0 



where A._, \ , ..., \ are real numbers, tt is the set of all affine 
combinations of X. X is said to generate tt. A k-flat in R is also 
known as a k-dimensional affine subspace of R . 
Equivalently , a k-flat tt is a set, 

tt = [x c R n | A x = b} (1-4) 

where A is a fixed real (n-k) x n matrix of rank n-k, and b is a 
fixed real (n-k)-vector. Both x and b are column vectors, and A x 
is the matrix product of A and x. (1-3) and (1-4) are equivalent 
under the assumption that Ax = b, for i =0, 1, ..., k. 

An (n-l)-flat is commonly known as a hyperplane . From (1-4) we see 
that a k-flat is an intersection of n-k independent hyperplanes. If k = n, 
the condition A x = b in (1-4) vanishes. Thus the only n-flat is R 
itself. A 0-flat is just a single real point. 



3 
Example 1-4 In R , 0,1,2,3-f lats are respectively: a single real 



3 
point, a line, a plane, and R itself. 



A 0-flat is a bounded set. A k-flat, where 1 < k ^ n, is an 
unbounded set. 



□ 
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Given k + 1 independent real points in R , there is one and only 
one k-flat containing these real points, and none of smaller dimension. 



1.1.4 Half-flats 



A half k-flat T in R , where 1 <; k ^ n, is a k-dimensional unbounded 
convex subset of R , 



T = {x g R n | A x = b and a x <. p} (1-5) 



where A x = b defines a k-flat, a is a fixed real n-vector (row) 
which is not in the row space of A, and p is a fixed real number. 
The set 

{x g R n | A x = b and a x = p} 
is a (k-l)-flat known as the boundary (k-l)-flat of the half k-flat t. 

3 

Example 1-5 In R , half 1,2,3-flats are respectively: a half-line, 

a half-plane, and a half-space. The corresponding boundary 0,1,2-flats 
are respectively: a point, a line, and a plane. 



□ 



A half n-flat, 



[x g R n | a x <; p} (1-6) 
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where n is a fixed real n-vector (row), cc ^ 0, and p is a fixed real 
number, is commonly known as a half-space . Its boundary (n-l)-flat 
is a boundary hyperplane. An intersection of a half -space (1-6) and 
a half-space, 

(x?R n | a ^P) (1-7) 

is the boundary hyperplane common to both (1-6) and (1-7). 

Let X be a subset of R . A half-space in R which contains X, and 
whose boundary hyperplane contains one or more real points in X, is 
known as a supporting half-space for X. Its boundary hyperplane 
is known as a supporting hyperplane for X. 



1.1.5 Convex Polyhedra 

A convex polyhedron P in R is a non-empty intersection of a 
finite number of half-spaces in R , 

P = (x g R n j A x <; b} 

where A is a fixed real m x n matrix and b is a fixed real m-vector 
(column). We shall say that P is an (n,k) -polyhedron if the rank 
of A is k. 

P may be bounded or unbounded, and may have any dimension between 
and n inclusive. Flats and half-flats are special cases of convex 
polyhedra. 
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The structure of a convex polyhedron is best understood by 
considering its faces. A face of a convex polyhedron P is a 
non-empty intersection of P with one or more supporting hyperplanes. 

Let a x <, b be any one inequality in the system A x ^ b defining 
P. The half-space 

T={ X fR n Jax <; b} 

may or may not be a supporting half-space for P. t is a supporting 
half-space if 

[x € P ! a x = b} (1-8) 

is non-empty. In this case, (1-8) is a face of P. 

Let S be an arbitrary subset of {1, 2, ..., m}. Let A denote 
the submatrix of A comprised of those rows of A whose indices appear 
in S, and let b denote the subvector of b comprised of those components 
of b whose indices appear in S. A subset P g of P, 

P s = {x g P | A s x = b s } 

if non-empty, is a face of P. Furthermore, every face of P may be 
obtained in this manner. 

If P g is a face of P, and we let S 1 , where S' c S, identity a maximal 

set of linearly independent rows of A , then P = P . Thus we may 

o S S 

generate all faces P g of P by considering only those sets of indices S 
which identify linearly independent sets of rows of A. We refer to such 
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sets of indices as face subsets . We include the case S = zp , where co 
is the empty set, so that P, which is equal to P , is a face of itself. 
Since there are only a finite number of ways to form S, P has a 
finite number of faces. The faces of P are themselves convex polyhedra. 

If a face of P has dimension k, we refer to that face as a k-f ace 
of P. 0, 1-faces of P, if they exist, are commonly known as vertices 
and edges of P, respectively. 

If P and P , arc faces of P and P c P , , then P is a subface 
of P ,. Clearly, if S' c S then P is a subface of P . . 

b b S 

A minimal face of P is one which has no proper subface. 
Let P be a convex (n ,k) -polyhedron , where < k < n. Hoffman 
and Kruskal [15] show that a face P of P is minimal iff A consists 

b O 

of k linearly independent rows of A. Furthermore, they show that if 
P is a minimal face of P, then 

p s = '' x f rI1 ! A s x = b s^ 

That is, all minimal faces of P are (n-k)-flats. 

The set of all faces of a convex (n ,k) -polyhedron P is partially 
ordered by c. The greatest member of this partially ordered set is P 
itself. There is no least member, in general. All faces fall into 
ranks, each rank being composed of all faces of a given dimension. 
The minimal faces of P, which are all (n-k)-flats, constitute the 
lowest rank. The next rank consists of all (n-k+l)-f aces of P, and 
so on. The highest rank consists of P alone, whose dimension lies 
between n-k and n inclusive. 
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Example 1-6 Fig. 1-3 shows graphical representations of several 

2 
convex polyhedra in R and partial ordering diagrams displaying 

the face structure of each. 



□ 



The boundary of a face of a convex polyhedron P is the union 
of all proper subfaces of that face. A boundary point of a face 
is a real point in the boundary of that face. A minimal face of 
P has no boundary; that is, its boundary is empty. The boundary 
of P itself is the union of all proper faces of P. 





3' p 




r {3? 



Ft.l Pw 



?M Pl« Pt37 F\»i Pu* Pt31 





Pji.*l Pn.** Pu.31 Pt«,2j Pia,u 



(a) 



(b) 



(c) 



(d) 



(e) 



Fig. 1-3 



20 

If P is a bounded convex polyhedron, then every face of P is 
bounded. Since every minimal face of P is flat, and the only 
bounded flat is a 0-flat, it follows that every minimal face 
of P is a 0-flat (or vertex), and that P is an (n,n) -polyhedron. 

A bounded convex polyhedron P is equal to the convex hull 
of its vertices (see Hadley [10]). Thus the vertices of P 
completely specify P. 

There is a similar result which applies to any convex polyhedron, 
bounded or unbounded. In order to state it, we first generalize 
the notion of a minimal face. 

A face of a convex polyhedron P may or may not be equal to the 
convex hull of its boundary. A face of P is said to be reducible 
if it is equal to the convex hull of its boundary, and irreducible 
if it is not. Flats and half-flats are irreducible faces. 

Klee [19] shows that a face of a convex polyhedron P is 
irreducible iff that face is a flat or a half-flat. He shows 
farther that P is equal to the convex hull of its irreducible faces. 

A minimal face of P is irreducible, but an irreducible face 
of P is not necessarily minimal. The converse does not hold 
because a half-flat is irreducible but not minimal. The notions 
of an irreducible face and a minimal face coincide in the case of 
bounded convex polyhedra since no face of a bounded convex polyhedron 
can be an (unbounded) half-flat. 

If an irreducible face of P is a half-flat, then its boundary 
is a flat which is also an irreducible face of P. A maximal irreducible face 
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of P is one which is not contained in another irreducible face of P. 
Klee [19] proves the following sharper result: 

A convex polyhedron P is equal to the convex hull of its maximal 
irreducible faces. 

Example 1-7 The minimal, irreducible, and maximal irreducible faces 
of the convex polyhedra represented in Fig. 1-3 are tabulated below: 



Fig. 


Bounded/ 


Minimal 


Irreducible 




1-3 


Unbounded 


Faces 


Faces 




(a) 


Bounded 


P 


P 




(b) 


Unbounded 


P {3] 


P > P (3} 




(c) 


Unbounded 


P fD 5 P {2} 


P {1}' P (2} 


P ! 


(d) 


Bounded 


P {1,2}> P {1,3}' 
P {2,3} 


P {1,2}' P {1,3}' 
P {2,3} 


p c 


(e) 


Unbounded 


P {1,2}' P (2,3) 


P (l}' P {3}' 
P {1,2}' P {2,3} 


P ( 



Maximal 

Irreducible 

Faces 

P 

P 



{1}' P {2] 



{1,2}' P {1,3} 



{2,3} 



P P 

{!}' {3} 



D 
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A class of convex polyhedra of special interest are those 
arising in connection with linear programming problems. A linear 
programming problem has the form: 

Maximize c x 

Subject to A x <; b (1-9) 

x ;> (1-10) 

where c is a fixed real n-vector (row), A is a fixed real m x n matrix, 
and b is a fixed real m-vector. The linear function ex is known as 
an objective function. A linear programming problem may instead 
involve minimizing an objective function. By simply negating c we 
may change such a minimization problem into a maximization problem. 
The set P of real points satisfying the inequalities (1-9) and (1-10), 

P = (x € R n | A x ^ b and x ^ 0} (1-11) 

is a convex polyhedron. We shall refer to P as a linear programming 
polyhedron. The condition x s gives P some special characteristics 
which we now examine. 

Suppose the inequalities (1-9) and (1-10) are consolidated and 
(1-11) is rewritten as 

P = {x f R n | A' x <; b'} 
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where A' contains as a submatrix the n x n identity matrix (negated). 
The rank of A' is n, and it follows that P is an (n,n)-polyhedron, 
all of whose minimal faces are 0-flats. 

What about the irreducible faces of P? It can be shown that 
{x g R | x ^ 0], known as the non-negative orthant of R , contains 
no 1-flat. It follows that the non-negative orthant of R contains 
no k-flat, for 1 ^ k <; n, and no half k-flat, for 2 ^ k ^ n. Since 
P is wholly contained in the non-negative orthant of R , P contains 
no k-flat, for 1 < k < n, and no half k-flat for 2 ^ k <, n. Thus 
the irreducible faces of P, necessarily being flats or half-flats, 
must be 0-flats or half 1-flats. 

P is equal to the convex hull of a finite number of 0-flats and/or 
half 1-flats. 



1.1.6 Integral Convex Sets 



Let J be the ring of integers. Let J be the n-dimensional 
module over J. We refer to x c J as an integer point or an 
integer n-vector . 

For any subset X c R , we define I(X) to be the set of integer 



points contained in X. 



I(X) = [x (z J n | x f_ X} 
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A convex subset X c R is said to be integral (or to have the 
integral property) iff X is the convex hull of the set of integer 
points contained in X, that is, iff X = H(I(X)). Since I(X) s X 
it follows that H(I(X)) c H(X) = X. Thus X is integral iff 
X c H(I(X)), that is, iff every x £ X is a convex combination of 
integer points in X. 

We are primarily interested in convex polyhedra having the 
integral property or, more simply, integral convex polyhedra. 
In Chapter 2 we develop necessary and sufficient conditions for 
a convex polyhedron to have the integral property. 

As the following example shows, not every integral convex 
set is an integral convex polyhedron. 
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Example 1-8 Let X = (x , x , x , . . . ) be the infinite subset 

2 
of J defined inductively as follows: 

x° = (0, 0) 

x n+1 =x n + (n+1, 1). 



Fig. 1-4 shows the first few members of X and a portion of the boundary 
of H(X). H(X) is an integral convex set because every x f H(X) is a convex 
combination of integer points in H(X). Clearly H(X) is equal to the 
intersection of an infinite number of half-planes and therefore is not 
a convex polyhedron. 

D 
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Fig. 1-4 



If X is a convex subset of R then there exists a unique convex 
subset X' of R n defined by X 1 = H(I(X)). In words, X 1 is the convex 
hull of the set of integer points contained in X. We point out two 
extreme cases. First, if X is an integral convex set then 
X' = H(I(X)) = X. Second, if X contains no integer points then 
X' = H(cp) = cp, where cp is the empty set. 

The following lemma is a direct consequence of our definitions. 

Lemma 1-1 

Let X and X* be convex subsets of R n . Then X 1 = H(I(X)) iff 



and 



(i) X 1 is integral 



(ii) I(X') = I(X). 
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Example 1-9 Consider the following system of linear inequalities 

(a) 4x <; 15 

(b) 2x <; 9 

(c) -3x £ -1 

The shaded interval of Fig. 1-5 (i) represents the set of real 
points x which satisfy (1-12). 



(1-12) 



Cc) 



Co.) a CW)- 



O l Z 3 4- 5 



C '•") 



O l Z 3 4- 5 



(i) 



(ii) 



Fig. 1-5 



Some of these real points are integer points and these are indicated 
with dots. In Fig. 1-5 (ii) the shaded interval is the convex hull 
of these integer points. It represents the set of real points 
which satisfy the following system of linear inequalities: 
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(a) x * |^J =3 



(b) -x < \=±\ = -1 



LtJ 



where I \J denotes the greatest integer less than or equal to 
the real number \. The shaded interval of Fig. 1-5 (i) is an 
example of a convex polyhedron. The shaded interval of Fig. 1-5 (ii) 
is an example of an integral convex polyhedron. The process of 
transforming the former into the latter is an example of integralization. 

□ 

Example 1-10 Consider the following system of linear inequalities: 



(a) 2x ± + 2x 2 <; 15 



(b) -4x x - 4x 2 <; -18 



(1-13) 



The shaded region of Fig. 1-6 (i) represents the set of real points 
x = (x.. , x ) which satisfy (1-13). It contains a countably infinite 
number of integer points, some of which are indicated with dots. In 
Fig. 1-6 (ii) the shaded region is the convex hull of these integer points. 
It represents the set of real points which satisfy the following system 
of linear inequalities : 

(a) x l + x 2 £ I — I = 7 

(b) - X]L - x 2 <; |j^-J = -5 

In Fig. 1-6 (i) we have an example of an unbounded convex polyhedron, 
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and in Fig. 1-6 (ii) an example of an unbounded integral convex polyhedron 
Again, the transformation of the former into the latter is an example 
of integralization. 



□ 
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Fig. 1-6 



Example 1-11 Consider the system 



(a) x 1 + 7x 2 <; 29 



(b) 4x x - 5x 2 <; 12 



(1-14) 



(c) -7x 1 - 4x 2 «j -27 
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Fig. 1-7 (i) illustrates the convex polyhedron defined by (1-14), 
while Fig. l-17(ii) illustrates the integral convex polyhedron 
which is the convex hull of integer points contained in the former, 





CO 



W) 



CO 



(i) 



(ii) 



Fig. 1-7 



The integral convex polyhedron of Fig. 1-7 (ii) is defined by the 
system of inequalities, 



(a) 



x„ < 



(b) x 1 -x 2 <_ 

(c) -x.. -x 2 <, 

(d) -x, < 
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If we were to imagine the x..-x„ plane as a peg-board, with pegs 
placed in all integer point positions, then a rubber band in tension, 
made to conform with the boundary of the shaded region of Fig. l-7(i), 
would, when released, assume the form of the boundary of the shaded 
region of Fig. 1-7 (ii), assuming that the rubber band remains in 
tension at the completion of its contraction. In its final state, 
the rubber band is supported entirely by pegs. Every vertex of 
the region it encloses is an integer point. This is a physical 
method for integralizing the bounded convex polyhedron of this 
example. 

In three dimensional space, a rubber sphere made to conform 
with the boundary of a bounded convex polyhedron and enclosing 
small marbles somehow fixed in integer point positions would do 
the job. In higher dimensional space it is difficult to imagine 
any physical process for accomplishing integralization. What is 
nedded, of course, is a mathematical procedure for performing 
the integralization transformation. 



In the chapters to come we shall be primarily concerned with 
two questions : 

(i) Given a convex polyhedron, how does one tell if it is an 
integral polyhedron? 

(ii) Given a non-integral convex polyhedron, how does one 
integralize it? 



□ 
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We answer the first question in Chapter 2 and we answer the 
second question for a restricted class of polyhedra in Chapters 
4, 5 and 6. 



1«3 Integer Linear Programming and Related Research 

In this section we seek to accomplish two goals. Our first 
goal is to establish the relevance of integral convex polyhedra 
and the process of integralization to integer linear programming. 
Our second goal is to describe some of the research done by others 
concerning integral convex polyhedra and integralization. 



1.3.1 Integer Linear Programming 

As we stated in the preceding section, a linear programming 
problem has the following form: 



Maximize 



c x 



Subject to A x ^ b > 
x s 



(1-15) 
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where c is a fixed real n-vector, A is a fixed real m y n matrix, 
and b is a fixed real m-vector. The inequalities in (1-15) define 
a convex polyhedron P. 

A commonly used algorithm for solving (1-15) is the simplex 
method developed by Dantzig [5]. The fundamental idea of the 
method is: if an optimal solution to (1-15) exists then some 
vertex of P is an optimal solution. Beginning with any vertex 
of P, the simplex method generates a sequence of vertices of P 
for which the objective function ex takes on non-decreasing values. 
The algorithm terminates when an optimum vertex is found, or 
gives an indication that no optimal solution exists. 

The simplex method has proven to be an efficient algorithm 
and has found many applications in operations research, management 
science and economics. 

An integer linear programming problem has the form: 



Maximize 



c x 



Subject to A x ^ b 
x ;> 
x an integer vector 



(1-16) 



where c is a fixed integer n-vector, A is a fixed m x n matrix, 
and b is a fixed integer m-vector. There are variations to (1-16), 
notably the more general formulation in which all constants are 
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allowed to be real. In another variation, only a fixed subset of 
the components of x are required to be integers; the remaining 
components are allowed to be real. We shall take (1-16) to be 
the standard integer linear programming problem. 

Again, the inequality constraints of (1-16) define a convex 
polyhedron P. The simplex method, when applied to (1-16), does 
not work in general simply because it is no longer true that if 
an optimal solution to (1-16) exists then some vertex of P is an 
optimal solution. Some or all of the vertices of P may not be 
integer points. Said in another way , optimal solutions to (1-16), 
if they exist, may be embedded in the interior of P and not 
accessible to any vertex generating algorithm. 

If P happens to be an integral convex polyhedron then, as we 
shall see in Chapter 2, every vertex of P is an integer point and 
the simplex method will produce an optimal solution to (1-16) 
for any objective function ex. 

Suppose P is not an integral convex polyhedron, as is usually 
the case. If we integralize P to obtain the integral convex 
polyhedron P', where P' = H(I(P)), then the simplex method applied to 
P' will produce an optimum solution to (1-16) for any objective 
function ex. We acknowledge the fact that there is a certain amount 
of naivety in this approach to integer linear programming. This 
stems from the fact that P' is potentially a very complex object — 
complex in terms of the number of inequalities needed to define P' — 



35 

and it certainly seems inefficient to construct P 1 in order to 
discover just one of its vertices. 

However, we offer the following observation, made by Gomory, 
in support of the idea of integralization as an important subject of 
research. We quote from Gomory [11] : 

"Since any algorithm for the integer programming problem, whether 
related to linear programming, branch and bound, exhaustive search, 
or whatever, must end up finding a vertex of P', information on 
P' seems relevant to any approach to the integer programming problem. 
Yet information about P' is very difficult to obtain." 
Although the last sentence of this quotation does not contribute 
to our argument, we include it to avoid mis-representation by omission. 

It has been our goal to gain some insight into the structural 
relationship between P 1 and P by focusing on the integral property 
and the process of integralization for convex polyhedra in general, 
not just linear programming polyhedra. We believe that such insight 
can contribute to the improvement of existing integer linear 
programming algorithms and the development of new ones. 



1.3.2 Related Research 

In this subsection we describe some of the known results related 
to integral convex polyhedra and integralization. 
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The major results concerning integral convex polyhedra known 
to us are contained in the paper by Hoffman and Kruskal [15]. 
They consider convex polyhedra defined by system of linear inequalities 
Ax ^ b in which all constants are integers, and define such a convex 
polyhedron P to be integral iff every minimal face of P contains 
an integer point. If the minimal faces of P are vertices then 
this definition says that P is integral iff every vertex of P 
is an integer point. As we shall see in Chapter 2, their definition 
is equivalent to our own. 

Hoffman and Kruskal obtain the following result concerning 
the integral property for convex polyhedra. Let A be an integer 
m y n matrix of rank k. Let S denote a subset of row indices 
of A and let A c denote the submatrix comprised of rows of A identified 
by S. Let gcd [A„] denote the greatest common divisor of all 
|S|-order minors of A„. Then Ax ^ b defines an integral convex 
polyhedron P(b) for every integer m-vector b iff 

gcd [A s ] =1 (1-17) 

for all k-subsets S of linearly independent rows of A. We observe 
that for a particular integer m-vector b, condition (1-17) is 
sufficient for P(b) to be integral but not necessary. 

Hoffman and Kruskal also obtain a result concerning the integral 
property for linear programming polyhedra. An integer m x n matrix 
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is said to have the unimodular property iff every minor (of every 
order) of A has value 1, 0, or -1. Entries in a unimodular matrix 
are therefore restricted to be 1, 0, or -1. They show that the 
inequalities Ax < b and x 2; define an integral convex polyhedron 
Q(b) for every integer m-vector b iff A is a unimodular matrix. 
We observe that for a particular integer m-vector b, the condition 
that A be unimodular is sufficient for Q(b) to be integral but 
not necessary. Unimodular matrices occur in the integer linear 
programming formulations of several combinatorial optimization 
problems, including the network flow problem and the maximal 
matching problem for bipartite graphs. 

Concerning integralization, Gomory in [10] are [11] gives 
a method for integralizing a corner polyhedron P — that is, 
a polyhedron defined by Ax ^ b where A is an integer n x n matrix 
of rank n and b is an integer n-vector. Basically the method 
involves generating a large but finite number k of auxiliary 
inequalities , 

ex 2: 1, i = 1, 2, . . . , k 

where c X , for i=l, 2, ...,k, is obtained as an integer solution 
to an equation over a finite Abelian group derived from the matrix A. 
Gomory shows that these inequalities define a polyhedron whose vertices 
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correspond one-to-one to inequalities defining P 1 , where P' = H(I(P)), 
and indicates how to obtain the latter from the former. Thus the 
problem of integralizing a corner polyhedron becomes one of 
finding all vertices of an auxiliary polyhedron, which can be 
accomplished by linear programming methods. 

We mentioned previously that the end product P' of integralizing 
a convex polyhedron P is potentially a very complex object. Rubin [22] 
obtains an interesting quantitative result concerning this complexity. 
Let a^, a„> a_, •••, a, , ... be the infinite sequence of integers 
defined by the recurrence, 



k = Vl + \-2 for k * 3 



a, = 1, a„ = 1. 



This is the well known Fibonacci sequence. Rubin defines an infinite 
sequence P^ P 2 , Pg, . . . , P fc , . .. of bounded (2,2)-polyhedra based 
on the Fibonacci sequence as follows. P, , for k i 1, is defined by 
the three inequalities , 

, 2 

a 2k X l + ^k+l X 2 * a 2k+l " i 



x^ ^0, x 2 s 0. 



p k > for k a 1, has 3 vertices and 3 edges. He shows that P' where 

P k = H ( I ( p k ))> for k a- 1, has exactly k+3 vertices and k+3 edges. Thus 
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for a 2-dimensional bounded polyhedron P, it may require 
arbitrarily many inequalities to define P' = H(I(P)), even 
though P is defined by just three inequalities. 
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CHAPTER 2 



THE INTEGRAL PROPERTY IN CONVEX POLYHEDRA 



In this chapter we develop several results concerning convex 
polyhedra having the integral property. Using our definition of 
an integral convex set X as one which is the convex hull of the 
set I(X) of integer points in X, we derive necessary and sufficient 
conditions for convex polyhedra to have the integral property. 

The chapter is divided into three sections. In the first 
we deal with integral flats, in the second with integral half- 
flats, and in the third with integral convex polyhedra in general. 



2.1 Integral Flats 

A k-flat tt in R is generated by taking all affine combinations 
of an independent set of k+1 real points. Lemma 2-1 tells us when 
tt is integral in terms of its generating set. 

Lemma 2-1 

A k-flat tt is integral iff rr is generated by an independent set 
[x , x , . . . , x j of integer points. 
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Proof 

We first prove necessity of the condition. Assume tt is an integral 
k-flat. Since w is a k-flat, tt is generated by an independent set Y of 
k+1 real points. Since tt is integral, each real point in Y is a 
convex combination of integer points in I(tt). The dimension of I(tt) 
is at most k. If I(tt) has dimension less than k, then Y has dimension 
less than k, which is a contradiction. Thus I(tt) has dimension k and 
tt contains an independent set of k+1 integer points which generates tt. 

In order to prove sufficiency of the condition, assume tt is 

1k 
generated by an independent set {x , x , . . . , x } of integer points. 

We shall show that every x f. tt is a convex combination of integer 

points in tt. It will then follow that tt is integral. 

Let x be an arbitrary real point in tt, 

1 k 

x = A- n x + \, x + • • • + K x 

where \ n , ^, , . ... K are real numbers and 
1 k 

\ n + \. + • • • + K =1. 
1 k 

Equivalently we may write 

x = x + k, (x -x ) + X (x -x ) + • • ■ + K (x -x ) . (2-1; 



Let functions h„: R -» J and h.. : R -» J be defined as follows. 
Let h n (>0 be the greatest integer less than or equal to \ and 
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let h-CX.) be the smallest integer greater than X... (Thus 
h 1 (\ i ) - h Q (\ ) = 1, for all real \ .) let {0,l} k be the set 
of all binary k-tuples, a representative member of which is 
denoted by b = (b.., b„, . .., b,). We claim that x is a 
convex combination of the 2 integer points, 



b 
x 



= x° + ^ b^ (\ t ) (x^x ), all b e {0,l} k 



i=l X 



in tt. (Geometrically, these 2 integer points are vertices of a 
k-dimensional cube in ^ and we claim that x lies in this cube.) 

In order to exhibit this convex combination we define functions 
g Q : R -» R and g- : R -♦ R as follows: 



g <V - ^(x.)^. 



We then define function g, : R -> R as follows 



»b ^ -ir^.V' 



i-i x 



We claim that 



I g^) x b (2-2) 

b e f°»i} 1 



x = 
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where 

£ g b (M £ 1 for all b f {0,l} k (2-3) 

and 

£ ^(N) = 1. (2-4) 

b £ {0,l} k 

In (2-2), \ = (\ , \_, ..., >>.,) is the real k-tuple of coefficients 
appearing in (2-1). 

In order to prove this claim we first observe that (2-3) holds 
since < g Q (?0 ■$. 1 and ^ g.. PO < 1 for all real numbers *... 

Next we prove (2-4) by induction on k. For k = 1 we have 



g (\ 1 ) + g 1 (\ 1 ) = h 1 (x 1 )-h Q (\ 1 ) = 1. 

Assume (2-4) is true for some k, 1 ^ k < n. Then 

k+1 

I [I \<\>] 

b <= {0,l} k+1 i_1 

Y k 

■ [ g o<w + *i<w] • b e f 0)1}k [ TT \^1 



= *0<>Sm) + 8 1<W 



= h i<W - h o<W = x - 
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Finally we prove that (2-2) is valid. We have 



^ A ^ b 

/ g b ( M x 



b p (0,1 



,k 



k 



V" 



L , g b (X) [ x ° + I h b ( V ^^] 
b p {0,l} k i= l x 



V 



b c !0,1) 



g. (\)1 x° + 
b 



V 



b c {0,1" 



gb (M 



k 

L-i b . 1 ! 

i=l X 



c° + V r v 



1=1 b . ro,i} k 



g b (X) h b ( vl (xi - x °)- 

i 



The expression within brackets may be expanded as 



V 



b c {0,1} 



k 



g h (M h (k.) 
b b t i 



k 

<L , L h 

b c {0,l] k J=l 



b . j _ b . l 
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k 
b g {0,l} k J =1 



g (X.) h (A..) + gl (X.) hl (\.) 



(h 1 (\ i )-\ i ) h (X.) + (X..h (\.)) h 1 (\ 1 ) 



= ^.(^(X.) - h (\.)) 



= X. 

L 



Thus 

b 6 (0,1)' 



1 

,k i=l 



= x 



and we have expressed an arbitrary x f n as a convex combination of 
integer points in rr- It follows that ^ is integral. 
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Let A be an integer (n-k) x m matrix of rank n-k and let b 
be an integer n-vector. Our next objective is to determine what 
conditions must be imposed on A and b in order for Ax = b to 
define an integral k-flat. The preceding lemma allows us to 
rephrase the question as: What conditions must be imposed 
on A and b in order for the solution space of Ax = b to contain 
an independent set of k+1 integer points. 
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We first introduce the following notation. If B is an integer 

r x s matrix of rank r, where 1 < r ^ s, then we define gcd [B] 

s th 
to be the greatest common divisor of all ( r ) r -order minors 

of B, where (£) is a binomial coefficient. Notice that if r = s 
then B is a square matrix and gcd [B] = |det Bj , while if r = 1 
then B is a row vector and gcd [B] is the greatest common divisor 
of all components of B. 



Theorem 2-1 

Let A be an integer (n-k) x n matrix of rank n-k and let b 
be an integer (n-k) -vector. Then 

A x = b 
defines an integral k-flat tt iff 

gcd [A] = gcd [A ib]. 



Proof 

Assume that A x= b defines an integral k-flat rr« By Lemma 2-1, 
■n is generated by an independent set of k+1 integer points. Let 

x be an integer point in tt. ^ en A x = b « l&t s denote an 

S 
arbitrary (n-k-l)-subset of column indices of A and let A be the 

(n-k) x (n-k-1) submatrix of A comprised of those columns of A 

identified by S. Since 
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det [A S i b] = det [A S j Ax°] 



S. 



k th 



it follows that det [A ib] is an integer combination of (n-k) -order 
minors of A. Thus 

det [A ! b]. 



Since this holds for any S, we conclude that 

gcd [A] = gcd [A i b]. 

Now we assume that 

gcd [A] = gcd [A ■ b] 

and prove that A x = b defines an integral k-f lat rr. 

According to Jacobson [18], there exists an integer (n-k) x (n-k) 
matrix R with det R = + 1 and an integer n x n matrix G with det C = ± 1 
such that 

A = R A C 
where A is an integer (n-k) x n diagonal matrix, 



A = 



*n-k 



Furthermore, Jacobson [18] proves that 



gcd [A] = gcd [Al 



It is easily seen that 



gcd [A] = a, a„ • • • a , . 
1 z n-k 



We may write Ax = b as 



R A C x = b 



or as 



A v = b 



where 



y = C x 
b = R _1 b, 



If we let C' be the following integer (n+1) x ( n +l) matrix 



C' 



+ -■ 
I 1 



we see that dot C ' =detC=+l and 



[A ; b] = R [A ; bl C 1 



49 



Again, it follows from Jacobson [18] that 



gcd [A i b] = gcd [A 1 b] 



We now have 



gcd [A] = gcd [A ; b ] 



?cd [A | gcd [A j bl 



from which we conclude that 



ecd 



gcd 



n-k 



n-k 



Thus 



a, a r 

1 j 



n-k 



°1 a 2 



a l a 2 



a l b 2 



n-k 



a a ,, • • • a , 

I z. n-k 



a l 3 2 



n-k 



and 



a i ! b i 



a 2 | b 2 



a n-k b n-k . 
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The k-f lat rr defined by A y = b contains the following independent 
set of k+1 integer points, 







V a i 




b 2 /a 2 




• 




• 




b , /a , 
n-k n-k 














• 




• 




• 


















V*i 




b 2 /a 2 




• 




• 




• 




- - 




b n-k /a n-k 




1 









• 




• 




• 









.. . 









E A 




b 2 /a 2 




• 




• 
• 




- - 




b . /a . 
n-k n-k 


• •• 







1 




. 




• 




• 














V a i 



V a 2 



b J* 1 

n-k n-k 






Therefore jt is generated by these k+1 integer points and, by 
Lemma 2-1, ^ is integral. Since the k-f lat tt defined by A x = b 
is the image of jj under the integer transformation c" , ^ is 
also an integral k-f lat, which is what we set out to prove. 



□ 



Corollary 2-1 

A x = b defines an integral k-f lat iff A x = b has an 
integer solution . 
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Proof 

The necessity of the condition follows from Lemma 2-1. 
To prove sufficiency, assume A x = b has an integer solution 
x . Then, as in the necessity part of the proof of Theorem 2-1. 

gcd [A] = gcd [A ib]. 

Thus, by Theorem 2-1, A x = b defines an integral k-flat. 
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Suppose that we have a system A x = b of n-k linear equations 
which defines an integral k-flat. The following corollary shows 
that any (n-k')-subset , k £ k' < n, of these linear equations 
defines an integral k'-flat. Let S be an arbitrary (n-k')-subset 
of {1,2,..., n-k}, where k ^ k' < n. Let A denote the submatrix 
of A comprised of those rows of A whose indices appear in S and 
let b c denote the subvector of b comprised of those components 
of b whose indices appear in S. 

Corollary 2-2 

If A x = b defines an integral k-flat then A Q x = b s defines 
an integral k'-flat. 



Proof 

Suppose A x = b defines an integral k-flat. By Corollary 2-1. 


A x = b has an integer solution x . But x is also an integer 
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solution to A x = b . Thus, again by Corollary 2-1, A x = b 



defines an integral k'-flat. 



□ 



2.2 Integral Half-flats 



We now consider the conditions under which a half-flat is 
integral. Let t be a half k-flat in R , 



t = {x(= r |Ax=b and a x ^ p} 

and let rr be its boundary (k-l)-flat. Again, we consider all 
constants to be integers. 

Theorem 2-2 

A half k-flat T is integral iff its boundary (k-l)-flat jj is 
integral. 

Proof 

Assume t is integral. Since t is not empty, I(t) is not empty. 
If rr contains an integer point then, by Corollary 2-1, ^ is integral. 
Otherwise 

a x < p for all x € I(t) 
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or 



a x <; p-1 for all x f- I( T ). 



Thus A x = b and a x <. p-1 define a new half k-flat t' such that 
t' c t a nd I(t') = I(t)« But this contradicts our assumption 
that t is integral. Thus rr contains an integer point and tt is 
integral. 

Now assume that rr is integral. We shall argue that T is 
the convex hull of a countably infinite number of integral 
(k-l)-flats, and thus t is itself integral. 

Let 

A 



A = gcd 



a 



and let rr. be the (k-l)-flat defined by 



a 



P-JA 



It can be shown that t is the convex hull of the set of (k-l)-flats 
{tt » tti» tt 2 » ••*}• We sha- 11 show that each (k-l)-flat ^ in the 



set is integral. 
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First of all, -r n is integral by assumption. We then have from 
Theorem 2-1 that 



gcd 



gcd 



A I b 



-X I P 



We next observe that each (n-k+1) -order minor of 



nr . P-JA 



which includes the right hand column differs from the corresponding 
minor of 



A 



b 



r - • 

a , P 



by a multiple of /\. It follows that 



gcd 



gcd 



A b 

*. I P~J A 
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and by Theorem 2-1, jj. is integral. Then, since each x £ T is a 
convex combination of real points in it. , tt, , TT„ , ... » each of 
which is a convex combination of integer points in T , it follows 
that each x f r is a convex combination of integer points in t« 
Thus T is integral. r-i 



2.3 Integral Convex Polvhedra 

We are now prepared to derive some conditions under which a 
convex polyhedra has the integral property. Let P be a convex 
polyhedron in R , 

P = {x e R n | A x £ b}. 

As before, we assume all elements of A and b to be integers. 

Theorem 2-3 

P is integral iff every face of P is integral. 

Proof 

Suppose every face of P is integral. Since P is the convex hull 
of its irreducible faces all of which are integral, P is itself integral. 
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Conversely, assume P is integral. If P is a flat then its 
only face is P itself which is integral by assumption. Otherwise 
P has one or more proper faces. Let P g be any proper face of P, 

p s = (xfP l A s x= V 

where S specifies a subset of the inequalities defining P. 

Let x g P Q , Since x £ P, x is a convex combination of integer 

12 m B 
points x , x , ..., x £ r, 



m m 

I c. x 1 , 0<c.<l, I c. -1. 
1=1 1=1 



We claim that x , x , ..., x f P , that is 



Otherwise 



A x 1 = b s , for i = 1,2,..., m. 



Ax < b„ for some i 



and 



m m 



A g x - A g 



1 c i xl - I c i A s xX < b s 



1=1 1=1 



which contradicts the fact that x g P g . Thus x is a convex combination 
of integer points in P g , and P g is integral. j-. 
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The following corollary gives a necessary and sufficient 
condition on the system of inequalities A x 5; b defining P for 
P to be integral. 

Corollary 2-3 

P is integral iff 



gcd [A g ] = gcd [A g ib g ] 



for all face subsets S. 



Proof 



Assume P is integral. By Theorem 2-3 every face P g of P is 
integral. Let jto be tne flat defined by A. x = b . Since P„ Q tt s > 
tt must contain an integer point, and by Corollary 2-1 and Theorem 2-1, 

gcd [Ag] = gcd [Ag i bg]. 

Now assume that 

gcd [Ag] = gcd [A g ■ b g ] 

for all face subsets S of P. Consider the irreducible faces of P. 
An irreducible face P_ of P is either a flat or a half-flat. If P g 
is a flat then P q is integral by Theorem 2-1. If P_ is a half-flat 
then its boundary flat is an irreducible face of P which is integral 
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by Theorem 2-1. Then, by Theorem 2-2, P is integral. Since P 
is the convex hull of its irreducible faces, all of which are 
integral, P is integral. i— -, 

We may further restrict the condition of the preceding corollary. 
This is done in the following. 

Corollary 2-4 

P is integral iff 

gcd [Ag] = gcd [Ag Jbg] 

for all minimal-face subsets S. 



Proof 

Necessity follows directly from Corollary 2-3. In order to prove 
that the condition is sufficient, assume that 

gcd [Ag] = gcd [A g ;bg] 

for all minimal-face subsets S. Consider the irreducible faces of P. 
An irreducible face P c of P is either: (i) a flat which is a minimal 
face of P, or (ii) a half-flat whose boundary flat is a minimal face 
of P. In the first case, P c is integral by Theorem 2-1. In the second 
case, the boundary flat of P is integral by Theorem 2-1, and thus P 
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is integral by Theorem 2-2. Since P is the convex hull of its 
irreducible faces, all of which are integral, P is integral. 



D 



Corollary 2-4 tells us that P is integral if all of its 
minimal faces are integral. If P is an (n,k) -polyhedron then 
the rank of A in A x ^ b is k and the minimal faces of P are all 
(n-k)-flats. If k = n then the minimal faces of P are all 
vertices. For example, this is the case for all linear programming 
polyhedra since the non-negative condition x £ insures a rank n 
coefficient matrix. In such cases, Corollary 2-4 asserts that P 
is integral iff every vertex of P is an integer point. 
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CHAPTER 3 



THE TOTALLY INTEGRAL PROPERTY FOR CONVEX POLYHEDRA 



In this chapter we define a new property which we call the totally 
integral property. We show that the class of totally integral polyhedra 
is a proper subclass of the class of integral polyhedra. We then 
consider the class of linear programming polyhedra having the totally 
integral property and show that this class can be characterized 
algebraically using some of the group theoretic ideas developed by 
Gomory [9 ] . 



3.1 Totally Integral Convex Polyhedra 

We begin by defining the totally integral property. Let P be 
a convex polyhedron defined by the system of inequalities A x ^ b, 
where, as usual, all elements of A and b are integers. Let S be 
a subset of row indices identifying a linearly independent set 
of rows of A. Then the system of equations A„ x = b defines a flat 
which we denote by tt s » We shall say that jt is a flat of P and that 
S is a flat subset Notice that a face subset is always a flat subset, 
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but that a flat subset need not be a face subset. It tt is a flat 

S 

of P then it is possible that P , where 



p s ■ p n tt s 



is empty and thus not a face of P. 

A convex polyhedron P is said to be totally integral iff 
every flat tt„ of P is integral. 



Example 3-1 Fig. 3-1 (a) illustrates a totally integral convex 
polyhedron. Fig. 3-l(b) illustrates an integral convex polyhedron 
which is not totally integral. 



□ 
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Fig. 3-1 
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The following lemma gives us a necessary and sufficient 
condition for P to be totally integral, expressed in terms of 
A and b. 

Lemma 3-1 

P is totally integral iff 



gcd [A_] = gcd [A |b ] 



S ' S" 
for all flat subsets S. 

Proof 

The lemma follows directly from the definition of the totally 
integral property and Theorem 2-1. 



□ 



Suppose matrix A has rank n-k, where ^ k < n. We shall say that 
a flat Tr ° f p is minimal iff S identifies n-k linearly independent 
rows of A. Thus a minimal flat of P is a k-flat. We may restrict 
the condition of the preceding lemma as follows. 



Lemma 3-2 

P is totally integral iff 



gcd [A s ] = gcd [A g jbg] 



for all minimal-flat subsets S. 
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Proof 

Necessity follows directly from Lemma 3-1. To prove sufficiency, 
assume 

gcd [A g ] = gcd [A g |bg] 

for all minimal-flat subsets S. Let R be any flat subset. We may 
augment R with additional rows of A, if necessary, to obtain a 
minimal-flat subset S such that R c S. By assumption, 

gcd [Ag] = gcd [Ag |b g ] 

and by Corollary 2-2 and Theorem 2-1, 

gcd [A^ = gcd [^ ;b R ]. 
It follows from Lemma 3-1 that P is totally integral. __ 

We see immediately that a totally integral convex polyhedron P 
is always integral. Every minimal face of P is a minimal flat of 
P. If P is totally integral then every minimal face of P is integral 
and, by Corollary 2-4, P is integral. 

The converse is not true — an integral convex polyhedron P 
need not be totally integral. It is possible that every minimal 
face of P is integral, but that one or more minimal flats of P 
which are not faces do not have the integral property. 
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Suppose we have an integer m y n matrix A and we let T(A) be 
the set of integer m-vectors b such that A x <; b defines a totally 
integral convex polyhderon. We show that T(A) forms a group under 
vector addition. 

We first observe that T(A) always contains the all zero m-vector. 
Applying Lemma 3-1 with b = 0, since 

gcd [Ag] = gcd [Ag ;0] 

for all flat subsets S, it follows that g T(A). Furthermore, T(A) 
is closed under vector addition. For if b and b are two integer 
m-vectors in T(A) then 

gcd [Ag] = gcd [A s ;bg] 
and 

gCd [A g ] = gCd tAgibg] 

for all flat subsets S, and it follows that 

gcd [A s ] = gcd [Ag | b* + bg] 

1 2 
for all flat subsets S. Thus b +b q T(A). Moreover, T(A) contains 

inverses for each of its members. If b f T(A) then 

gcd [Ag] = gcd [Ag ; b g ] 

and 

gcd [Ag] = gcd [A s i -bg] 
for all flat subsets S. Thus -b f T(A). 
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In the following section we completely characterize the group 
T(A) for linear programming polyhedra. 



3.2 Totally Integral Linear Programming Polyhedra 

We now focus our attention on those convex polyhedra contained in 
the non-negative orthant of R — that is, convex polyhedra P defined 
by 

A x ^ b 

(3-1) 
x s; 

where again A is an integer m y n matrix and b is an integer m-vector. 
Given matrix A, let T(A) denote the set of integer m-vectors b for 
which (3-1) defines a totally integral convex polyhedron. In the 
following, we obtain a complete characterization of T(A), using 
group theoretic methods originally developed by Gomory [9]. 
Specifically, we show that T(A) may be characterized implicitly as 
the set of integer solutions to a single equation expressed over 
a finite Abelian group. 
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Writing the constraints (3-1) defining P as 



A 






b 





X 


< 


— 


I 










where I is an n x n identity matrix, we see that the coefficient matrix 
has rank n. It follows that all minimal flats of P are 0-flats 



(or points), as are all minimal faces of P. From Lemma 3-2, P is 
totally integral iff all minimal flats of P are integer points. 
Also, from Corollary 2-4, P is integral iff all minimal faces of P 
are integer points. 

We reformulate the constraints (3-1) defining P in the following 
equivalent form, 



[Ai I] 



x 



= b 



(3-2) 



X £ 0, Z 2; 



(3-3) 



where I is an m y ni identity matrix and z is an m-vector of slack 
variables. A basic solution to (3-2) is one which is described 
as follows. Let B be any basis of [Ail], i.e. a non-singular 
m y m submatrix of [A i I], and let x be the corresponding 



m-subvector of 



Let N be the remaining m x n submatrix 



of [A i I] and let x be the corresponding n-subvector of 
We may write (3-2) as, 



67 



The unique solution 



[B ; N] 



~B 

*N 



A B 



= b. 



to (3-4) given by 



(3-4) 



x B -B b 



*N 



= 



is a basic solution to (3-2). If x s then the basic solution 

-i B 

X, 



B 



is a basic feasible solution to (3-2), and B is a feasible 

basis for (3-2). Notice that whether or not a basis B of [Ail] is 
feasible depends upon the particular right hand side vector b in 
(3-2). 

It is easily seen that 0-flats of P correspond to basic solutions 
of (3-2) and that 0-faces (or vertices) of P correspond to basic 
feasible solutions to (3-2). Thus, P is totally integral iff every 
basic solution to (3-2) is an integer point, and P is integral iff 
every basic feasible solution to (3-2) is an integer point. 

We now take a slightly different point of view in expressing 
the condition for P to be totally integral. Let B be a basis 
and let M(B) denote the set of all integer combinations of columns 



of B. Then the basic solution 



B 




is an integer point iff the 
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solution x to Bx = b is an integer point or, equivalently, iff 
B " 

b £ M(B). There are at most ( n + m ) bases of [A | I] , which 
we denote by B^ B 2 , ..., B,. It follows that P is totally 
integral iff 

b f M(B jL ) for i = 1,2, ..., k 

or, equivalently, iff 

k 
bg f) M(B .). 
i=l 

Recalling our defining of T(A) as the set of integer m-vectors 

b for which P — defined by (3-2) and (3-3) with A fixed — • 

is totally integral, we see that 

k 

T(A) = f| M(B .). 

1=1 

We note that M(B i ) , for i = 1,2, ..., k, forms an Abelian group 
under vector addition. Moreover, this group is a subgroup of 
the Abelian group M(I) under vector addition, where I is the 
m x m identity matrix. It follows that T(A), as an intersection 
of subgroups of M(I) , is also a subgroup of M(I). We summarize 
these results in the following. 
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Lemma 3-3 

Given a fixed integer m x n matrix A, the set T(A) of all 
integer m-vectors b for which the constraints 



[A il] 



= b 



x ;> 0, z £ 

define a totally integral convex polyhedron is a subgroup of M(I). 
Moreover 

k 
T(A) = fl M(B.) 
i=l 

where B. , for i = 1,2, ..., k, are all the bases of [A ; I]. _ 

Before proceeding with our characterization of T(A), we should 
point out that a similar development is not possible for the integral 
property. The fundamental difficulty with the integral property 
is its connection with feasible bases. Which bases of [Ail] are 
feasible depends on the particular choice of b as a right hand side 
vector. Given an integer m x n matrix A, suppose we let Q(A) be 

the set of integer m-vectors b for which (3-2) and (3-3) define 

1 2 
an integral convex polyhedron. Then, if b and b are in Q(A), 

1 2 
it is not true in general that b + b is in Q(A). This is due 
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to the fact that the feasible bases associated with b may not be 

2 
the same as those associated with b . Thus Q(A) is not closed 

under vector addition. 

The real advantage of dealing with the totally integral property, 
and our motivation in defining it, is that it involves all bases 
of [A ; I] for every choice of right hand side vector b. As a result 
it leads to an interesting algebraic characterization of T(A). 

We begin by considering a single basis B of [A i I] and showing 
how M(B) is characterized as the set of integer solutions to an 
equation expressed over a finite Abelian group whose structure 
is completely known. An integer m-vector b is in M(B) iff B b 
is an integer m-vector. If we let 

A = |det B| 

then B b is an integer m-vector iff 

AB~ b - mod A (3-5) 

where the columns of AB are integer m-vectors, is the all zero 
m-vector, and the congruence applies componentwise. Consider the 
group M(AB ). If we introduce the congruence relation on M(AB ) 
which relates two integer m-vectors iff they are congruent modulo A 

componentwise, then the quotient group with respect to this congruence 

, . M^'B" 1 ) 
relation is , where I is the m v m identity matrix. Moreover, 

M(AE) 
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if we let p. , p_, ..., P denote respectively the congruence 
1 z m 

classes containing each of the columns of ^B , and if we 
now let denote the congruence class containing the all 
zero m-vector, then (3-5) may be expressed as the following 

mCab" 1 ) 

equation over the group G = ' . . _. , 

M.\al) 

m 



£ b. p. = 0. (3-6) 



i=l 

Thus an integer m-vector b is in M(B) iff b satisfies the group 
equation (3-6) over G. 

We next show that (3-6) can be replaced by an equation over 
a new group which is isomorphic to G, whose size and structure 
is more easily revealed. First we note that if H is any 
group isomorphic to G, with isomorphism cp: G -t H, then the 
equation (3-6) over G is satisfied by some integer m-vector 
b iff the following equation over H 

m 



£ * ± cp(P £ ) = cp<°) 



1=1 
is satisfied by b. 

m( a b" ) 

We claim that the group G = — is isomorphic to the 

M( A I) 

group G 1 = i.ft\ • The non-singular linear transformation -B is 
an isomorphism from M( A B ) to M(I) and from M(^I) to M(B). It 
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follows that — B takes congruence classes in G into congruence 
A 

classes in G 1 and that this mapping from G to G' is an isomorphism. 

We now take the isomorphism one step further. Smith [25] has 
shown that by performing certain elementary row and column operations 
on matrix B one obtains a diagonal matrix B, 



B = 







■m 



where f 1 , £„, . .., c are positive integers, €,"?2* *"* ^m ~ ^ > 
and f.|£.,, f°r i = 1,2, ..., m-1. Matrix B may be expressed 
in terms of matrix B as follows , 



B = R B C 



where R-and C are integer m x m matrices with Idet R| = |det C| = 1. 

We show that the group G' = : i , is isomorphic to the group 

G" = -c" . Since Idet Rl = 1, R is an isomorphism from M(I) to M(I). 
M(B) ' 

Moreover, R is an isomorphism from M(B) to M(RB). But M(RB) = M(B) , 
which can be shown as follows. Let y g M(B). Then y = RBCx = RB(Cx). 

for some integer m-vector x. Thus y £ M(RB) and M(B) c M(RB). Next 

-1 - -1 
let y f M(RB). Then y = RBx = RBC(C x) = B(C x) , for some integer 

m-vector x. Since |det Cl = 1, C x is an integer m-vector and 

therefore y f M(B). Thus M(RB) c M(B) , and M(RB) = M(B). It follows 
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that R is an isomorphism from M(B) to M(B). Furthermore it follows 

that R takes congruence classes in G 1 into congruence classes in G", 

and this mapping from G 1 to G" is an isomorphism. 

The structure of G" = \~( is apparent. If we let J, denote 

M(B) v k 

the group of integers {0, 1, 2, ..., k-1} under addition modulo 

k, then it is readily seen that G" is isomorphic to a direct product 

group H, where 



H = J x J _ X • • • X J • 
fe l ^2 %i 



The isomorphism from G" to H maps a congruence class in G" into the 

integer m-tuple obtained by taking any representative member of 

the class and reducing its i component modulo £., for i = 1, 2, ..., m. 

The order of H is £..•£„• ... •£ = A* Although we have expressed 

H as an m-fold direct product, in many cases several of the £. 's have 

value 1, in which case the corresponding components of the direct 

product contribute nothing to the structure of H and can be 

omitted. 

The following diagram summarizes the sequence of isomorphic 
groups we have been discussing. 



1 reduction 

A B R modulo C^, f 2 , ..., fO 



m(ab _1 ) 



Mill ~ Hill ~ T T T 

^^ = ^^ = J x J X • • • x J 



M(aD M(B) m(B) <=! X (= 2 x '•• x ^ 
G G 1 G" H 
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We began with the equation (3-6) over the group G as a means 
of characterizing the members of M(B) . We may now replace that 
equation with an equation over the group H. Letting cd denote 
the isomorphism from G to H, we replace equation (3-6) with 



m 

) b. h. 
l l 

i=l 



(3-7) 



where 



h. = o((B.) 



- RB(Ab ). mod 



R. mod 
l 



e 2 



th 
(The subscript i on a matrix denotes the i column of that matrix.) 

Equation (3-7) is the desired characterization of M(B). An integer 

m-vector b is in M(B) iff b satisfies equation (3-7). 

Now we are ready to characterize T(A), the group of all m-vectors 

b for which the constraints (3-2) and (3-3), with matrix A fixed, 

define a totally integral convex polyhedron. We recall from Lemma 3-3 

that 



T(A) = Pi M(B.) 
i=l L 
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where B. , for i = 1, 2, ..., k, are all the bases of [Ail]. 

For each basis B. we have seen that M^) may be characterized 
as the set of integer solutions to an equation 

m 

Y b . h . . = 
L J ij 

over a finite Abelian group E^ of order ^ = |det B ± \. It follows 

that T(A) may be characterized as the set of integer solutions to 

k simultaneous equations, 

m 
> b. h. = 0, for i = 1, 2, ..., k 

j=l 
over the finite Abelian groups E^, H 2 , ..., E^. Equivalently, 
T(A) may be characterized as the set of integer solutions to a 
single equation, 

m 

Y b. k. = 

j-l 
over a finite Abelian group K, where 

K = H x x H 2 x ' • • X \ 

and 

kj = (h lj9 h 2j h kj ), for j = 1, 2, ..., m. 

The order of K.is A-i*A 2 * ••• "A^« 
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We summarize these results in the following theorem. 

Theorem 3-1 

Given a fixed integer m x n matrix A, the set T(A) of all 
integer m-vectors b for which the constraints 



[All] 



x 



= b 



x ^ 0, z ;> 

define a totally integral convex polyhedron is equal to the set of 
integer m-vectors b which satisfy an equation 

m 

Y b. k. = 
L J J 

over a finite Abelian group K derived from all bases B. , B , ..., B 
of [A i I]. The order of K is A 'A-* ... *A , where 



A . = |det B.|. 



□ 



One conclusion we may draw from Theorem 3-1 is that T(A) never 
consists of just the all zero m-vector. If we let o(k.) denote the 
order of k in the group K, and if we let be the diagonal m x ™ matrix, 
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°(k 1 ) 



o(k 2 ) 



o(k ) 
m 



we see that 



M(0) c T(A). 



Tlius T(A) contains a countably Infinite number of m-vectors. 

In terms of convex polyhedra having the totally integral property, 
these results are somewhat surprising. They tell us that for any 
matrix A there exist a countably infinite number of right hand 
side vectors b which give rise to totally integral polyhedra. 
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CHAPTER 4 



TRANSFORMATION OF (n ,k) -POLYHEDRA INTO (k ,k)-POLYHEDRA 



The general problem of Integralizing convex (n,k) -polyhedra 
may be approached inductively with respect to one of the parameters 
n or k. Recall that an (n ,k)-polyhedron is defined by Ax <, b , where 
A is an m x n matrix of rank k, 1 <, k <; n. In this chapter we show 
that the class of (n,k)-polyhedra for fixed k and arbitrary n, 
n ;> k, have a great deal in common with respect to integralization. 
In particular we show that an (n,k)-polyhedron P , where n 2: k, may 

be transformed into a (k,k)-polyhedron P by a linear transformation 

n k 
T: R _» R which preserves integralization. 

In order to explain what we mean when we say that T preserves 

integralization, we refer to the diagram below. If we let ' denote 



K 
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the operation of integralization, then this diagram expresses the 
fact that 



T(P') = T(P )'. 
^ n n 

Thus in order to integralize P to obtain P', we may instead integralize 

n n 

P, to obtain P' and then use the inverse transformation T to obtain 
k k 

P'. (A word of explanation concerning T is in order. As we shall see. 
n 

T is a many-to-one transformation. Thus it does not make sense to talk 

-Ik n 

about T (y), y f R , as being a single real point x € R • Rather, by 

T (y) we shall mean 

T _1 (y) = {x <= R n | T(x) = y}. 

With this interpretation of T we may say that T (P') = P'.) 

Thus the purpose of this chapter is to show that the general 
problem of integralizing (n,k)-polyhedra reduces (somewhat) to the 
problem of integralizing (k,k)-polyhedra. In the first section of 
this chapter we show how to obtain a suitable transformation T for 
a particular (n,k) -polyhedron P . In the second section we show 
that such a transformation preserves integralization. And in the 
third section we begin our inductive approach to integralization 
by examining the simple case, k = 1. 
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4.1 Existence of a Suitable Transformation 

Let P n be an ( n > k ) -polyhedron defined by Ax ^ b, where A is an 
integer m x n matrix of rank k, and b is an integer m-vector. Any 
integer kxn matrix T of rank k whose rows span the row space of A, 
when viewed as a linear transformation from R n to R , transforms 
P n into a ( k > k )~P ol y he dron P defined by Cy ^ b, where y = Tx and 
C is an m x k matrix (not necessarily integer) of rank k defined 
by the matrix equation A = CT. However, we require a matrix T 
having two additional properties : 

(i) gcd [T] = 1 

(ii) every row of A is an integer combination 
or rows of T. 

In the following section we show that a matrix T having these 
properties preserves integralization. In this section we are 
concerned with the existence of such a matrix T. 

The conditions we have imposed on T are dependent only on 
the matrix A and are independent of the vector b. The following 
lemma assures us that a matrix T exists having the required properties. 
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Lemma 4-1 



If A is an integer m x n matrix of rank k then there exists an 
integer k y n matrix T of rank k such that (i) gcd [T] = 1 and 
(ii) every row of A is an integer combination of rows of T. 



Proof 



Let U be any k x n submatrix of A having rank k. According 
to Jacobson [18], U may be written as 

U = R U C 

where R is an integer k x k matrix with det R = + 1, C is an integer 
n x n matrix with det C = + 1, and 

U = [U 1 ; 0] 

where U' is a diagonal k x k matrix and is the k x (n-k) zero matrix. 
Let T = [I i 0], where I is the k x k identity matrix and let 

T = R T C. (4-1) 

It is clear that gcd [f] = 1. Jacobson [18] shows that gcd [T] = gcd [T] = 1. 

We now show that every row of A is an integer combination of rows 
of T. Let a be any row of A, Since A has rank k and U is a k x n submatrix 
of A having rank k, the rows of U span the row space of A. The row a 
may be written uniquely as a linear combination of the rows of U, 
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a = cU 
where c is a rational k-vector (row). Then we may write 

a = cU = c(RUC) = cR(U'T)C = c(RU'r" 1 ) (RTC) 

= c'T 

where c 1 is a rational k vector (row) given by c ' = c(RU'R -1 ). 
We show that c 1 must be an integer k-vector. 

We may write c' = |_c'J + <c '> where J_c'J denotes 
the greatest integer function applied componentwise to c 1 and 
<c'> is the vector of fractional parts of c'. Each component 
of <c'> is a non-negative rational number less than 1. 
Now, since a = c'T = l_c ' J T + <c '> T and Lc ' J T is an 
integer n-vector, we see that <c '> T is also an integer n-vector. 
We shall deduce from this that <c '> = 0. 

We represent <c '> uniquely by 

<c > = - (e^ e 2 , ..., e^) (4-2) 

where d is a positive integer, e^ e 2 , ..., e, are non-negative 
integers, <j e,. < d, and gcd (e.^ e 2> ..., e fc , d) =1. Let S 

be a subset of k column indices of T such that the k x k submatrix 

S 
T comprised of the columns identified by S has non-zero 

determinant. Since <c *>T is an integer m-vector we see that 
g 
<c'>T is an integer k vector, say a*. Solving the system 
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g 
of equations a ' = <c '> T for <c '> we find that <c '> has the form 

— (f - , f , ..., f , ) , where f 1 , f _ , ..., f, are integers. 

|det T b | l l k 1 2 k 

Comparing this with (4-2) we conclude that d J det T . Similarly, 

we may show that d divides all k yk minors of T. Thus d I gcd [T] = 1, 

and d = 1. It follows that e n = e„ = . . . = e, =0 and <c '■> = 0. 

12 "k ^ □ 



4.2 The Transformation Preserves Integralization 

The transformation T whose existence is guaranteed by Lemma 
4-1 transforms the (n,k) -polyhedron P defined by Ax ^ b into the 
(k,k) -polyhedron P, defined by Cy ^ b, where y = Tx and A = CT. 
Property (ii) of T insures that C is an integer matrix. We shall 

write P k = T(P n ) and P n = t" (P ). In general, for a subset 

k -1 
Y e R , T (Y) will have the meaning 

T _1 (Y) = {x f R n | Tx g Y}. 

We argue next that T takes faces of P into faces of P, . 

n k 

Let (P ) be an arbitrary face of P , 
no n 

( Vs = {x * P n I A s x - V 
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and let (P,) be the face 



(p v ) s - (y f.K ! c q y - b s }. 



- k / s i; *- * k I ~s 



Since A = CT we see that A g = C T. It follows that T((P n ) g ) = (P k > s - 

It also follows that T" 1 ((P k ) s ) = (* ) g . Tbe minimal faces of P n 

are (n-k)-flats in R n . Their images under T are minimal faces of 

R . 0-flats in R . 
k' 

Let ' denote the operation of integralization. Suppose we 
have computed P' an integral (k,k)-polyhedron defined by C'y fib 1 , 
where C 1 is an integer m' y k matrix of rank k and b' is an integer 
m'-vector. Then T (P') is an (n,k) -polyhedron defined by A'x <; b', 
where A 1 is an integer m' x n matrix given by A* = C'T. The following 
lemma asserts that the transformation T preserves integralization. 



Lemma 4-2 



P* = T" 1 (P.') 
n k 



Proof 

We show that t" 1 (P^) is integral and I(t"" 1 (P ] ^)) = I(P n ). It will 

follow from Lemma 1-1 that T _ (P,' ) = P'. 

k n 



r " 1(Pi 



We first show that T (P') is integral. Every minimal face of 



T (P') is the image of a minimal face of P' under T and conversely, 
k K 

Minimal faces of P' are vertices. Since P,' is integral, its vertices 
k k 
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are integer points in R . Let y be any vertex of P'. Then 
T _1 (y ) is a minimal face of t" 1 (P^). T _1 (y ) is the (n-k)-flat 
defined by, 

T^Cy ) = {x g R n I Tx = y°}. 



Here we make use of property (i) of T, namely gcd [T] = 1. 
It follows that 

gcd [T] = gcd [T ; y ] 

and, by Theorem 2-1, T (y ) is an integral flat. It follows that 

every minimal face of T (P') is integral, and by Corollary 2-4, 

T (P') is integral. 

We now show that I (T~ (P')) = I(P ). Since P' c P. and 

k n k k 

T" 1 (P .) = P we have t" 1 (P 1 ') c P and I(t" 1 (P 1 ')) c I(P ). 
k n k n k n 

Consider the set difference I(P ) - I (T ( p ^))- Assume there 
is an integer point x in this set difference. Since x f I(P ), 
Tx = y is an integer point in P. , hence in P.'. But then 

J or ^> ^ 

x f. I(T (P,')) which contradicts our assumption. We conclude 

that I(P ) - I(T" 1 (P^)) is empty. Thus I(T~ 1 (P 1 |.)) = I(P n >. 

Having shown that T~ (P') is integral and that I(T ( p ^)) = 

I(P ), it follows from Lemma 1-1 that T~ (P.') = P'. _ 

v n ' k n [j 



Summarizing these results, we have seen that given a system 
of linear inequalities Ax <, b defining an (n,k)-polyhedron P , 
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there exists a transformation T which transforms P into the 
(k,k)-polyhedron P, defined by Cy <; b, where y = Tx and A = CT. 
Furthermore, if we integralize P, to obtain P' defined by 
C'y <; b', then T _1 (P') defined by A'x <; b 1 , where A 1 = C'T, 
is equal to P'. Thus, without loss of generality, we may 
restrict ourselves to the problem of integralizing 
(k,k)-polyhedra. 

We examine the case k=l in the following section to 
illustrate the transformation and integralization of 
(n-l)-polyhedra. 



4.3 Transformation and Integralization of (n,l)Polyhedra 

Let P be an (n, 1) -polyhedron defined by Ax <. b , where A 

n 

is an integer m x n matrix of rank 1, and b is an integer m-vector. 

Let T = — a., where a. is any non-zero row of A and A- = gcd [a.]. 
Aj J J J J 

Then gcd [T] = 1. Moreover, every row a. of A, i = 1, 2, ..., m, 
is an integer multiple of T, a. = c. T, where c. is an integer equal 
to + gcd [a.]. Thus T satisfies the conditions for an integralization 
preserving transformation. 

T(P ) is the (l,l)-polyhedron P defined by Cy ^ b , where 



n 



1 



y = Tx and C is the integer m x 1 matrix defined by A = CT. 

According to Lemma 4-2, T (P,') = P'. 

° In 
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The problem of integralizing P 1 is easily solved. Consider 
a single half-line T . defined by c^y <; b ± > c ± $ 0. We integralize 
T . as follows. Let Z. be the half-line defined by 



-l*J 



if c ± > 



■y < [4J if c * < °- 



i 



The boundary point of T . is an integer so that, by Theorem 2-2, 
t. is integral. Furthermore, I(t.) = K.T.)- It follows from 

X it 

Lemma 1-1 that t- = r[ = HClOi^)). 

Among all those inequalities c.y ^ b for which c > 0, if 

any, let the j inequality be one for which the ratio is 

i 
smallest. Similarly, among all those inequalities c y < b for 

which c. < 0, if any, let the k inequality be one for which 

b. 
the ratio — is smallest. Let P. be defined by the following 
" C i 1 

system of at least one and at most two inequalities. 



y < 



■y < 
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We see that I (P.) = I(P-). If 



b, 




b. 


k 


> 


JL 


~c, 




c . 


k 




J_ 



then 



I (Pi) = I(P-i) = 9, where cp denotes the empty set. Thus 

P^ = Pj = cp. Otherwise, by Corollary 2-4, P, is integral. 

It follows from Lemma 1-1 that P- = P ' . 

T (P') is an (n,l) -polyhedron defined by the following 

system of at least one and at most two inequalities, 

b. 

_-L 

c . 

J 

b. 



TX <; 



-Tx <; 



~c, 



By Lemma 4-2, t" 1 (P') = P'. 

In 



Example 4-1 Let P. be the (3,1) -polyhedron defined by 



-8 4 12 
12 -6 -18 
-4 2 6 



13 
7 
9 



Let T = ± [-8 4 12] = [-2 1 3]. Then P ± = T(P 3 ) is the 
(1,1) -polyhedron defined by, 



mm _ 




— *■ 


4 




13 


-6 


y <: 


7 


2 




9 


■. _ 




_ — 
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W 



e integralize P- and obtain the (1 ,l)-polyhedron Pj defined by. 



y < 



■1.T,. 



Then T (P ' ) is the (3 ,1) -polyhedron defined by 



*2 

2 



1 3 
-1 3 





— — 




x l 




x 2 




X 3 




— - 



and T _1 (Pp = P^. 



D 
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1 i 
■j- a -x s 



d. 

l 



for i = 1, 2. Since I( T !) = I(T.) we have 

' 1 1 



i(p) = i(t 1 n t 2 ) = i(t^ n Tp 



so that P may be replaced by t' q t'. 

In the following we assume that this first integralization step 
has been performed on P. That is, we assume that P is defined by 
Ax <. b where 



A = 



a ll a 12 



a 21 a 22 



, b - 



and 



gcd(a lis a 12 ) = 1 



gcd(a 21 , a 22 ) = 1. 



Let A = det A. We assume that the order of the two inequalities 
Ax < b is such that A > 0. 

Our goal is to generate the system of inequalities A'x ^ b' 
defining P 1 . We accomplish this by generating all normal vectors 
and vertices of I(P). In the following we define the notion of 
a normal vector of I(P) and a vertex of I(P). We then introduce 
a transformation of co-ordinates which makes the properties of 
normal vectors more accessible. In particular, it allows us to identity 
images of normal vectors as being atoms of a certain partially ordered set. 
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We then develop a two-dimensional generalization of the division 
theorem for integers and show how it is used to generate the 
finite set of these atoms. We then describe a second level 
two-dimensional division process involving these atoms which, 
together with the inverse transformation, generates all normal 
vectors and vertices of I(P). These lead directly to the system 
A'x^b 1 defining P'. 



5.1 Normal Vectors and Vertices 





that x is a boundary point of I(P) iff there exists a non-zero 

/cr 1N 



integer vector a = ( J such that 



£T«(x-x ) <; for all x £ I(P) 



where . denotes vector inner product. We shall say that CT 

is a support vector of I(P) at x . The inequality a. (x-x ) <; 

defines a supporting half-plane for I(P). 

We shall say that a boundary point of I(P) is a vertex of I(P) 

iff there exist two linearly independent support vectors of I(P) 


at x . 
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It is a straightforward matter to compute some of the boundary 

points of I(P). To do so we make use of some basic results in number 
theory. 

Let tt be an integral line defined by a«x = b, where a = ( ) 



is an integer vector, b is an integer, and gcd (a,, a„) = 1. 

. • , . 

An integer solution x to a*x = b may be determined as follows. 

The Euclidean algorithm is used to compute a (non-unique) integer 

vector r = ( r ) such that a«r = 1. Then x = br is an integer 

solution to a*x = b, and a-x = b may be written as a- (x-x ) = 0. 

The following lemma is stated without proof. (For a proof 
see, for example, Niven and Zuckerman [21].) For convenience 
in stating it, we let 



T = 



1 
■1 



Lemma 5-1 



If tt is a line defined by a* (x-x ) = where gcd (a.. , a ) = 1, 

then I(tt) is the set of integer points expressible as x + kTa, 

where k is an integer. 

□ 

Consider the integral boundary lines tt- and tt of the integral 
half-planes T^ and i^ defining P. tt, and tt 2 are defined by a 1 .x = b 
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2 i ' 

and a -x = b„ respectively, where a = ^ 



/3- • -l \ 

) and gcd (a..., a.„) = 1, 

.a.,,/ i-L i z 



for i = 1,2. The Euclidean algorithm gives two integer points x and 

x 2 in tt and tt 2 respectively. The defining equation for rr., may then 

be written as aMx-x 1 ) = 0, for i = 1,2. Lemma 5-1 tells us that 

I(tt ) is the set of integer points expressible as x + k Ta , where 
i *- 

k. is an integer. 
i 

2 
Certain integer points in Id^) satisfy a «x <, b 2 and thus are 

integer points in I(P). It is easy to show that KjT^) o I(P) is 

the set of integer points expressible as x -k^Ta , where k^ is 



an integer and 



k l * 



u 2 1 
b -a «x 



Similarly, I(tt ) n I(P) is the set of integer points expressible 

2 2 
as x + k„Ta , where k 2 is an integer and 



k 2 < 



b -a -x 



Both I(rr 1 ) n I(P) and I(tt 2 ) n I(P) are sets of boundary points of I(P). 
If x° is a boundary point of I(P) then we shall say that a non-zerc 



. 



integer vector CT is a normal vector of I(P) at x iff: 
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\ 
(i) a- (x-x°) < for all x c I(P) 



(ii) x° + To c I(P) ^ (5" 1 ) 



(iii) gcd (a , a 9 ) = 1. 



The inequality o,(x-x ) < defines a half-plane t which is a supporting 
half-plane for I (P) . The boundary line t defined by o (x-x ) = 
is a supporting line for I(P). u contains at least two integer 
points, x and either x + Tq or x - To", in I(P). The condition 
gcd (o\ , a 9 ) = 1 is included to provide a canonical form for a. 
If this condition were not imposed, then every positive integer 
scalar multiple kcr of a normal vector c would also be a normal 
vector, and as we shall see, this would be inconvenient for our 
purposes „ 

Two normal vectors of I(P) are evident. For every boundary 

point x in I(t,) ^ I(P), a is a normal vector of I(P) at x . 

2 2 

Similarly, for every boundary point x in I(tt 2 ) n I(P)» a is 

2 
a normal vector of I(P) at x . 

Since normal vectors are support vectors, we see that a boundary 
point x of I(P) is a vertex of I(P) if there exist two linearly 
independent normal vectors of I(P) at x . 

In the remainder of this chapter we develop a method for 
generating the system of inequalities defining H(I(P)). The method 
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involves generating all normal vectors and vertices of I(P). 
Clearly, if x is a vertex of I(P) and a is a normal vector 
of I(P) at x , then the inequality cr» (x-x ) ■£ defines a 
half-plane whose boundary line contains an edge of H(I(P)). 
It follows that the system of all inequalities of the form 
e«(x-x ) ^ 0, where x is a vertex of I(P) and a is a normal 
vector of I(P) at x , defines H(I(P)). 

We begin by introducing a transformation of co-ordinates 
which makes the properties of normal vectors more accessible. 



5.2 A Hull-Formation Preserving Transformation 



Let us regard the matrix A of the system Ax ^ b defining 

2 2 2 

P as a linear transformation A: R -> R . For all x g R , 

A: x -► y where y = Ax. Since A ^ 0, A is a non-singular 

linear transformation. Thus the inverse transformation A 

exists. A is given by 



-1 



1 
A 



a 22 " a 12 



l 21 



l ll 



A is not an integer matrix in general. 
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2 
We shall denote the image A(X) of a subset X c R under A by X 

so that 

» 2 
X = f y ^ R ! y = Ax, x c x} . 



Since 



P = <> c R [ Ax <: b} 



it follows that 



2 



P = [y c R I y < b}. 

Letting M(A) denote the set of all integer combinations of columns 
of A we see that 

I(P) X = {y a M(A) 1 y <; b}. 

We claim that the transformation A preserves the operation 
of forming the convex hull of I(P). This is illustrated in the 
following diagram. 



I(P) 



KP) 



H(I(P)) - H(I(P) ) 



This diagram expresses the fact that 

H(I(P)) X = H(I(P)*) 
which we prove in the following lemma. 

Lemma 5-2 

H(I(P))* = H(I(P)*) 

Proof 

We first show that H(I(P)) W c H(I(P)*). Let y c H (I(P))*. 

Then 



/ -* , 1 ? 

y = A( A. Q x + X x + X x ) 



where ,\ Q , X^ X,, > 0, X Q + ^ + ^ = 1, and x \ x \ x 2 e I(p) 

We then have 

1? 

y = \ Q Ax + A. Ax + X Ax 

0,1 2 

- * Q y + X^ + \ 2 y 

1 2 
where y , y , y ,: L(P) . Thus y c H(I(P) ). 

We next show that ii(I(P)") c H(I(P))*. Let y € H(I(P)*). 

Til en 

y - V° + \y l + \y 2 



where \ Q , X ± , ^ 2 (), ^ + ^ + ^ = ^ and y ^ ^ y 2 ^ i(p) ^ 



We then have 



99 



y = ^-qAx + X,Ax + ^ 2 Ax 



= A(a. q x + \.x + \_x ) 



A 1 O ^ 

where x , x , x £ I(P). Thus y £ H(I(P)) . 



D 



We should point out that unless A = ± 1» A does not preserve 

'k 'k 

the operation of integralization. That is, (P 1 ) £ (P ) * , where 

1 denotes the operation of integralization. In order to see this, 

* * 
consider P . P is a (2,2)-corner polyhedron defined by y ^ b. 

* * 

Since the vertex y = b of P is an integer point, P is integral . 

That is, P = (P )'. However, if A t ± 1, then it is not true 



in general that the vertex 



x = A d o 



f P is an integer point, 



in which case P is not integral. That is P' jt P. It follows 
that (P 1 )* # P = (P*)\ 



Example 5-1 Let P be the (2,2)-comer polyhedron defined by 



4 3 
-3 1 



23 
-5 



Fig. 5-l(a) represents portions of P, I(P) and H(I(P)). Fig. 5-l(b) 

□ 



represents portions of P , I(P) and H(I(P) ). 
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It can be shown that the linear transformation A preserves 

the face structure of H(I(P)). The image of a face of H(I(P)) 

* 

is a face of H(I(P) ) and conversely. 

Let x be a boundary point of I(P) and let a be a normal 
vector of I(P) at x . Then <j« (x-x ) ^ for all x £ I(P). In 
terms of I(P) we have 

o-'A'^y-y ) <; (5-2) 

for all y € I(P) , where y = Ax is a boundary point of I(P) . 
Multiplying (5-2) by A > to clear fractions we obtain 

a. (AA _1 ) (y-y°) <; 



or in terms of matrix multiplication only, 



c T (AA" 1 )(y-y°) <: o. 



Letting 



we have 



T T -1 
Y = cr (AA ) ( 5 _3) 



v(y-y ) £ 



for all y £ I(P) . We shall refer to v as a normal vector of I(P) 



at y . 

* 
Since normal vectors of I(P) play a key role in our development 

we shall examine their properties in some detail. First of all, 
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we see that a normal vector y of I(P)" is related to the corresponding 
normal vector a of T (P) by equation (5-3), or equivalently , taking 



the transpose of each side, by 



-1 T 
Y = (A A ) a . 



-IT 
Letting A = (AA ) we have 



Y = Ac 



(5-4) 



A is the cofactor matrix of A and is given by 



22 



-a 



21 



" a 12 a ll 



Letting M(A) denote the group of all integer combinations of columns 
of A under vector addition, we see from equation (5-4) that y g M(A) . 

The defining conditions (i), (ii) , and (iii) in (5-1) for a 
normal vector n of I(P) at x may be recast as defining conditions 
for a normal vector y of l(P)" at y as follows. If y is a boundary 



po 



int of I(P)" then a non-zero vector y 6 M(A) is a normal vector 







of I(P) at y iff: 



(i) y(y-y ) <, for all y G I(P)' 



(ii) y° ± Ty c I(P)' 



(5-5) 



(iii) X y i- M(A) for all real X, < X < 1. 
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We have already seen that condition (5-5) (i) follows from 
(5-1) (i). Condition (5-1) (ii) tells us that 



v* 

y + ATo- € I(P) . 



(5-6) 



Now AT = TA, as is easy to verify. Thus ATa = TAa = Tv and (5-6) 
becomes condition (5-5) (ii). Finally (5-1) (iii) implies that \a 
is not an integer vector for all real X, < "K < 1. Condition 
(5-5) (iii) follows. Thus, if CT is an integer vector which 
satisfies (5-1), then y = Aq is an integer vector which 
satisfies (5-5). The converse argument can also be made. 

Example 5-2 Let P be the (2,2)-convex polyhedron defined in 

Example 5-1. Vertices and associated normal vectors of I(P) 

* 
and I(P) are tabulated below. 



I(P) 


vertices 


© 


© 


(0 


normal 
vectors 


(?). (?) 


Cx), (!) 


©. G) 


* 

I(P) 


vertices 


(3) 


(2) 


(-U) 


normal 
vectors 


V13), uoj 


(£). © 


©. (^ 



104 



Notice that tabulated normal vectors of I(P) are images of the 

1 3 



corresponding tabulated normal vectors of I(P) under A = 



-3 



□ 



5.3 Normal Vectors as Atoms 

In this section we show that normal vectors of I(P) are 
necessarily atoms of a certain partially ordered subset of 
M(A). We begin with the following lemma. 

Lemma 5-3 

* 
If Y is a normal vector of I(P) at y then y 2; 0. 



Proof 



According to (5-5) (i) 



Y«(y-y°) <: for all y e I(P)*. 



Recall that I(P) is defined by 



(5-7) 



I(P) = {y € M(A) | y <: b}, 



Since y £ I(P) , y f M(A) and y ^ b. Consider y = y - ( n )« 

( Q J is an integer vector, rM f M(A). Thus y 1 £ M(A). 



Since A 



D 
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1 1 * 

Also, since y <. b and A > 0, y <: b. Therefore y € I(P) . 

Setting y = y in (5-7) we obtain -Y-, A <. or 

A similar argument using y = y - ( . j reveals that 

Y 2 * °- 

A normal vector y is, by definition, a member of M(A). 
If we let M(A) be defined by 

M(A) + = {z £ M(A) | z ;> 0} 

then Lemma 5-3 tells us that y is a member of M(A) . As we see next, 
a much stronger statement about y can be made. 

Consider the partially ordered set (M(A) , ^) . The ordering 

- + 
relation ^ is, as usual, taken componentwise. The element £ M(A) 

is the least element or universal lower bound. 

An atom of (M(A) , <;) is an element a € M(A) such that a £ 

and 



z^a =* z = or z 



for all z g M(A) + . 



Theorem 5-1 



•k 

If Y is a normal vector of I(P) at y then y is an atom of 



(M(A) + , <; ). 
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Proof 

Let y be a normal vector of I(P) at y . By definition, 
y £ M(A) and y j= 0. By Lemma 5-3 y £ M(A) + . 

Assume there is a z f_ M(A) , distinct from and -y such 
that z <; y. According to (5-5)(ii), y + Ty £ I(P)"„ We show 
that y + Tz and y + T(y-z) are also members of I(P) ' and use 
(5-5) (i) to conclude that y. (+Tz) = 0. It will follow that 

z = \y, for some real X, < \ < 1, which contradicts (5-5) (iii). 

* 

We first show that y + Tz f I(P) . Since z £ M(A) 

and TA = AT, it follows that Tz £ M(A). Since y £ M(A) , 

y + Tz e M(A). Now since y° + Ty c I(P)" we have y° + Ty <; b. 

Then since ^ z ^ y we have y + Tz < b . Thus y + Tz 6 I (P) \ 

* 

We next show that y + T(y-z) f I(P) . Since (y-z) f M(A) 

and TA = AT, it follows that T(y-z) f M(A). Since y <= M(A) we 

have y + T(y-z) f M(A) . Again y + Ty <; b. Since ^ (y-z) ^ y 

we have y ±T(y-z) ^ b. Thus y + T(y-z) f I(P)'\ 

Now we argue that y- (+Tz) = 0. Setting y = y + Tz in (5-5) (i) 
we obtain y. (+Tz) ^ 0. Setting y = y + T(y-z) in (5-5) (i) we 
obtain y. (+T(y-z)) < 0. Since y. (+Ty) = 0, this gives -y. (+Tz) <; 0. 
Thus y • (+Tz ) = . 

Finally y-(+Tz) = implies that z = Xy for some real number X. 
Furthermore O^z^y, z^0 and z ^ y, imply that < A. < 1. But 
z = ^Yj < A. < 1, contradicts (5-5) (iii). 



□ 
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ft 
Theorem 5-1 tells us that normal vectors of I(P) are necessarily 

atoms of (M(A) , ^). As we shall see, there are at most A + 1 such 

atoms. In the remaining sections we show how to generate these 

ft 
atoms and then show how they are used to compute H(I(P) ), and 

thus H(I(P)). 



5.4 A Two-Dimensional Division Theorem 

ft 4- 

Having shown that normal vectors of I(P) are atoms of (M(A) , ^), 
our next objective is to generate these atoms. To this end we present 
the following theorem, which may be viewed as a generalization of 
the division theorem for integers. It applies to integer 2-vectors. 

The notation ( f , j means a. > b- and a 2 < b«. All quantities 
are integers. 

Theorem 5-2 



Mi v /U, 



Given ( J and f j such that 

C 1 ) > C 2 ) 



> ,u„ v > 



/ U 3\ 
there exist unique q and ( J such that 
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U1 ) 


= q 


C 2 ) 




f 3 ' 


■V 


Vv 2 / 




S 



and 



U 2\ > /V * ° 



o < w 2 y < \v 3 



Proof 

f U 3\ 
We first show that q and I exist. 

Let q and u be determined by the division theorem applied to u^ 
and u„, 

^ = q u 2 - u 3 , u > u 3 ;> 0. 

Since u.. > u„ > 0, it follows that q > 1. 
Then v„ is given by 

v 3 = q v 2 " V l 
and since q > 1 and <; v, < v 2 , it follows that < v 2 < Vy 

f U 3\ 
We next show that q and ( J are unique. 

Assume 

("l) . , ("^ . ( U 3) (5-8) 

Vv.,/ M \v 2 / Vv / 
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and 



r 



^) . q < r>) . n) (5 . 9) 

1' X V Vv 3' 



where 



U 2\ > / u 3\ * ° / u 2\ > / U 3\ * ° 

) ) and . Z ) ( ■?) (5-10) 

< w 2 ■ < W / < \v 2 / < \v'y 



Subtracting (5-9) from (5-8) gives 



« - o fi 2 w:v dv (5-iD 



,v 2 / w 3 / w 3 



Now from (5-10), 



u 2 > u 3 > (5-12) 



and 



;> -u^ > -u 2 . (5-13) 



Adding (5-12) and (5-13) gives 



u 2 > u 3 - u 3 > -u 2 



which, together with the following equation from (5-11) , 



(q - q r ) u 2 = u 3 - u 3 ' 



implies 

q = q' 
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This, together with (5-11) implies 

CP " (vp . D 

/ U l\ f U 2\ 
Given two ordered pairs of integers ( j and ( J such that 

/ U l\ > / U 2\ > ° 
_ ( 1 ( J successive application of Theorem 5-2 yields 



'V ^ ^ v 2 
a sequence of ordered pairs 



,u > /^ > / U 3N > ... > /V2N > (VlN > (\\ 

o ^ Vvj^y < \v 2 ; < \v 3 ; < < w k _ 2 / < \v k-1 / < w k / 



= 



(5-14) 



The sequence (5-14) terminates when u, = 0, which occurs for some 

finite k since u- , u~ u, is a monotonically decreasing sequence 

and u- is a positive integer. 

There is some symmetry in the statement of Theorem 5-2 which 

allows us to apply the theorem in the reverse direction. If we 

/ u 2\ / u 3\ 
display the ordering relations among the components of I J and ( J 

as follows , 



\U, ' < \u_/ 



< -3 

we find that we have two ordered pairs to which Theorem 5-2 may be 
applied producing unique q and I J such that 



Ill 



and 



Q) = ^-C$ 



C 2 ) > C 1 ) * ° 

< Vu 2 ; < \a J 



Example 5-3 Let f J = V 3 ) ' (v ) = (s)* Successive application 
of Theorem 5-2 gives the sequence 

(?) 0) (A) (A) (i) Q- 

Applying Theorem 5-2 in the reverse direction, we may extend this 
sequence to the following one: 

(T) ft 1 ) (?) G) (A) (A) &) (X) . D 

The sequence of first components in (5-14) 

u l > u 2 > *** > \ = ° 

is one which would be produced by the Euclidean algorithm applied to 
u, and u_. The set of common divisors of u- and u„ is equal to 
the set of common divisors of u. and u . _ , for i = 1, 2, ..., k-1. 
In particular, if u. and u„ are relatively prime, then so are u. and 
u. ., for i = 1, 2, ..., k-1. 
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The sequence of second components in (5-14) 

V k> V k-l > '" > v l*° 

is one which would be produced by the Euclidean algorithm applied 
to v^ and v, - . Similar remarks about the set of common divisors 
of v jL and v,. +1 for i = 1, 2, . . . , k-1 apply here. 

Any three consecutive ordered pairs in the sequence (5-14) 
are related by an equation, 



0) - % c-3 - c::) 



i+r i+2' 

for i = 1, 2, . .., k-2. These equations may be used to express any 
ordered pair in the sequence (5-14) as an integer combination of 
any two consecutive ordered pairs in the sequence. 



Example 5-4 In the previous example, f „ ] may be expressed 
integer combination of (_, j and L_J as follows, 



as an 



(^ 



) - 7 Q - («)• 



Similarly, f , j may be expressed as an integer combination of 
( 1 J and ( J as follows , 

A) - "= (?) ♦ » ( 1 3 6 )" 

□ 
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5.5 Generation of Atoms 

In this section we show how to generate the set of atoms of 
(M(A) , <j). Before doing so, we establish some preliminary results 
concerning M(A). 

Recall that A is given by 



*22 



l 12 



*21 



*11 



where gcdCa^, a 12 ) = 1 and gcd(a 21> a 22 > = 1. For simplicity 
in our notation, we let 



a 



11 



a 



21 



a 12 
a 22 



where a xl = a 22 , a 12 = " a 2l» a 21 = " a 12 and a 22 ~ a ll* Thus 
gcd( ct , 11 , a 12 ) = 1 and 8 cd (a 21 » a 22 ) = *• 



Let a 1 = ( 11N ) and a 2 = ( 21 Y Th e Euclidean algorithm may 
\a 12 / Vx 22 / 

1 / r llN 2 / r 21\ 1 1 

be used to compute an r = ( ) and an r = ( J such that a «r = 1 

\r 12 / \r 22 / 



2 2 



and a «r = 1. They are not unique. Let g, be an integer such that 
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1 2 

g, s a -r mod a> i g. < i (5-15) 



and let g„ be an integer such that 



2 i 
g 2 = a «r mod A, «j g 2 < A (5-16) 



where A = det A = det A > 0. 

Integers g. and g_ are uniquely defined even though vectors 

12 2 

r and r are not. To see that g. is uniquely defined, let r 

2 2 

and (r )' be any two vector solutions x to the equation a *x = 1. 

2 2 2 

We may write a >x = 1 as ct • (x-r ) = 0. From Lemma 5-1, 

2 2 2 
(r ) ' = r + kToc for some integer k. It follows that 

12 12 
a • (r )' =a «r mod a» Thus g 1 is uniquely defined by (5-15). 

A similar argument reveals that g„ is uniquely defined by (5-16). 



Lemma 5-4 



g 1 g 2 s 1 mod A 



Proof 



gj^ g 2 = (a t ) (a «r ) mod £ 



= [V-r 1 ) (a 2 -r 2 ) - (cc^Ta 2 ) (r^Tr 2 )] mod A 
= [ 1 - A O^'Tr 2 )] mod A 
= 1 mod A 



a 
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1 2 
Now, let (r )' and (r )' be any two vector solutions to the 

1 2 

equations a «x = 1 and a *x = 1 respectively. Let 



where 



and let 



r 1 = (r 1 )' + k^a 1 



2 , L, 
a «(r )'-g. 



k l = 



2 2 2 
r = (r )' + k 2 Ta 



where 



1 , 2., 
g-,-a «(r )' 

k 2 ; — 



11 22 12 21 

Then a «r = 1, a «r = 1, a t =g 1 and a *r = g 2 - We define 

matrices C. and C„ as follows : 



C l = 



a 22 r 21 



" a 21 r 22 



C„ = 



r ll " a 12 



r i2 a ll 



We observe that det C- = 1 and det CL = 1. Let G- = AC 1 and 
G„ = AC 2 . Then 



G l = 



A Si 



1 



and G„ = 



1 



g 9 A 



(5-17) 
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Lemma 5-5 



M(A) = M(G ) = M(G 2 ) 



Proof 



We 



prove that M(A) = M(G 1 ) . Let z £ M(A) . Then z = Ax for 



some integer vector x. It follows that z = A(G n C, )x = AC, (C, x) , 



, x « 1 " 1 )x = AC 1 (C 1 " 1 
where C 1 is an integer matrix because det C. = 1, and C- x is an 



integer vector. Thus z £ M(AC, ) = M(G, ) and we have shown that 
M(A) c M(G 1 ). 

Now let z £ M(G..). Then z = G-x = (AC. )x = A(C.x) for some 
integer vector x. Since C..X is an integer vector, z £ M(A). 
Thus M(G 1 ) c M(A), and we have shown that M(A) = M^). 



A similar argument gives M(A) = M(Gn)« 



D 



Lemma 5-6 



gcd (A, g x ) = 1 and gcd (A, g 2 ) = 1. 



Proof 



We prove that gcd (A, g, ) = 1. From (5-17) we see that 



- 1 ) c M(G-). By Lemma 5-5, M(G ) = M(G ) , so that ( l ) £ M(G ). 
j = G x for some integer vector x. This says that 



s r 



1 = x 1 A + x 2 g 1 



from which we conclude that gcd (A, g-i)l !• It follows that 



?cd (A, g x ) = 1. 
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A similar proof shows that gcd (A, g ? ) = 1. 



We are now ready to generate the set of atoms of the partially 
ordered set (M(A) , ^). Suppose we apply our two-dimensional 
division theorem successively to the columns of G. . That is , 
suppose we let 

Q) - Q) - ft) 

and apply Theorem 5-2 successively to obtain the finite sequence 

o.c;):o:(5:-:(^):(t;):ct)- -> 

where ki2. We shall show that the set of integer vectors 
appearing in this sequence is precisely the set of atoms of 
(M(A)\ <). 

For convenience we let W be the 2 ^ k matrix, 



»- r x 



u, u„ u 



2 U 3 V-2 Vl Vi 
V l V 2 v 3 V k-2 Vl V 

/U 



/ \ 
whose columns w = ( J, i = 1, 2, ..., k, are the integer vectors 



appearing in the sequence (5-18). Thus 



w.ro, !" 3 ... V2 vi \i 

'- X ' v 3 v k-2 v k-l V 
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We first argue that every column w of W is a member of 
M(A) . As we noted in the previous section, the equations 
used to derive the columns of W allow us to express every 
column w as an integer combination of the columns of G- . 
Thus w £ M(G ). From Lemma 5-5, w g M(A). Furthermore, 



J, w £ m'" v+ 
We next show that 



since u. ^ and v. £ 0, w £ M(A) . 



Lv k-1 V L s 2 4j 



We begin by showing that u, , = 1. From Lemma 5-6, gcd (u.. , u_) = 1. 
Furthermore, as previously noted, gcd (u. , u .) = gcd (u- , u_), 
for i = 1, 2, ..., k-1. Thus gcd (vl . , u.) = 1. Since u^ = 
and u, _. > 0, it follows that u, . = 1. 

We now show that v, = A by expressing w as an integer 
combination of the columns of G- . The equations used to derive 
W are as follows : 

1 2 3 
w = q-w - w 

2 3 4 
w = q w - w 



k-3 k-2 k-1 
w = q, w - w 
k-3 

k-2 k-1 k 
w = \-2 W " w * 



(5-19) 



From the last of these equations we have 
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k k-1 k-2 
w = q „ w - w 



Since il = and ii - = 1, we see that \_ 2 = \_2' We then 
write (5-20) as 



k k-1 k-2 
\-2 W " "k-l*" 



Then using (5-21) and equations (5-19) we have 



k , k-2 k-3 N k-2 



w" = \-2 ( \-3 W " W } " Vl W 

, k-2 k-3 
= ( %-3\-2 " Vl )w " \-2 W 



From the last of these equations we see that 



v k = Ul v 2 - u 2 v x 



= A-l - g 1 '0 . 



Thus v, = A» 



(5-20) 



w = u, w " u, ,w . (5-21) 



k-2 k-3 

U k-3 W " \-2 W 



k i+1 i (5-22) 

w=u.w-u., n w v/ 

1 1+1 



= u„w - u w 



= u w - u„w . 
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k-1 



Finally we argue that v., = g„. Since w £ M(G.) we may 



x 1 A + gl x 2 



k-1 _ rt k-1 

write w = G.x for some integer vector x. Thus w 

1 \ x 

from which we conclude that 



Vi s Vk-i mod A - 



Multiplying both sides by g_ we obtain 



8 2 U k-l = g l 8 2 V k-l m ° d A * 
Using Lemma 5-4 and the fact that u, , = 1, we have 



V k-1 ~ g 2 mod A 
and this, together with «g v, , < Aj implies v , = g„. 
Thus we may write matrix W as 



u„ 



u. 



V 3 \-2 
W may also be generated from right to left by starting with the 

columns of G and successively applying Theorem 5-2 in the 

reverse direction. 



Example 5-5 Let 



103 -46 
-27 61 
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Then A = det A = (103) (61) - (-27) (-46) = 5041. The Euclidean 
algorithm gives 

21(103) + 47(-46) = 1 
9(-27) + 4(61) = 1 



Then 



and 



g x = 9(103) + 4 (-46) mod 5041 
= 743 mod 5041 



g 2 = 21(-27) + 47(61) mod 5041 
= 2300 mod 5041. 



Thus 



G l = 



5040 743 
1 



G„ = 



1 
2300 5041 







Matrix W is generated using either G or G and Theorem 5-2. 

5041 743 I 160 57 11 9 7 5 3 ' 1 
I I 

W = I I 

1 | 7 34 95 536 977 1418 1859 | 2300 5041 



We require two additional lemmas before we prove that the set 



□ 



of columns of W is the set of atoms of (M(A) , ^). 
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Lemma 5-7 



Let G. , i = 1, 2, 



k-1, be the 2 y 2 submatrix of W 



given by, 



G. = 



u. u. _ 
l l+l 



1 1+1 



Then M(G.) = M(A). Furthermore, det G. = det A = A. 
i ' l 



Proof 

The equations used to generate the columns of W allow us to 
express the columns of G. as integer combinations of the columns 
of G^. Thus M(G.) c M(Gj). The same equations allow us to 
express the columns of G- as integer combinations of the columns 
of G ± . Thus MCG^ c M(G jL ), and we have shown that M(G.) = MCG-) . 
By Lemma 5-5, M(G.) = M(A). 

Furthermore, from equation (5-22) we see that 



v, = u.v. ., - u. , ,v. 
k l i+I i+l i 



and since v, = A, we have det G. = A. 



D 



We define M (W) to be the set of all non-negative integer 
combinations of columns of W, 



m (w) = {z € J ! z = wx, x e J , X s: 0}. 
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Lemma 5~8 



M + (W) = M(A) + 



Proof 

Since each column of W is a member of M(A) , it follows that 
every non-negative integer combination of columns of W is a member 
of M(A) + . Thus M + (W) c M(A) + . 

We now argue that M(A) + c M + (W). Let z 6 M(A) + . If z = then 
clearly z € M (W). Therefore we can assume that z ^ 0, that is 
z. > or z. > or both. 

Consider the quantities 



Z l U i 



z 2 v. 



for i = 1, 2, ..., k. For i = 1 we have 



z l u l 



z 2 v x 



Z l A 



z 2 



= -z A <; 



(5-23) 



and for i = k we have 



l l \ 



Z 2 V k 



2 1 ° 



z 2 A 



= z 1 A i 0. 



(5-24) 



We claim that 
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z l 


u. 

1 


< 


Z 2 


v. 

1 





z l u i+l 



Z 2 V i+1 



(5-25) 



for i = 1, 2, ..., k-1. To prove this we observe that 



£ v ± < v i+1 and u ± > \i ±+1 ;> 0. 



(5-25) follows from this and our assumption that z. > or z„ > or both. 
Now, from (5-23), (5-24) and (5-25) we conclude that either: 

(i) for some i, 1 < i < k 



z, u. 
1 l 



Z 2 V i 



< and 



Z l U i+1 



Z 2 V i+1 



> 



or 



(ii) for some i, 1 ^ i ^ k 



Z l U i 



z v. 
2 i 



= 



We consider case (i). By Lemma 5-7, M(G.) = M(A) . Since z £ M(A) 
it follows that z £ M(G.) and z = G.x for some integer vector x. We show 
that x- > and x > 0. From Lemma 5-7, det G. = A- Applying Cramer's 
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rule, we solve z = G.x for x: 

1 



1 

x i = Z" 



U i+1 



'i+1 



> 



X 2 A 



u. 

l 



V. 



1 

A 



> 



Thus z is a non-negative integer combination of columns of G. . 

In case (ii) , similar reasoning gives z = x- ( j where x. 

i 
is a positive integer. 

+ - + + 

Thus we have shown that z £ M (W) . It follows that M(A) c M (W) . 



We conclude that M + (W) = M(A) + . 



D 



We are now ready to prove the main result of this section. 



Theorem 5-3 



The set of columns of W is the set of atoms of (M(A) , ^) . 



Proof 



We first show that every column w of W is an atom of (M(A) , ^) 



i + i - + 

We know that w £ M(A) and w ^ 0. Assume that there is a z £ M(A) 



126 

distinct from and w such that z ^ w . By Lemma 5-8, z may be 

expressed as 

k 

2 = I C i wl 
i-1 

where c. is a non-negative integer, for i = 1, 2, ..., k. If c^ = 0, 
for i = 1, 2, ..., k, then z = contrary to assumption. Thus 
c. > for some j, 1 ^ j g k, and it follows that vr ^ z. Since 
w J £ z and z ^ w 1 it follows that w 3 s w 1 , which implies that 
w^ = w , because 

U l > U 2 > *•* >\ 

v l < V 2 < •-< V 

Now we have w^ ^ z < w and w = w . We conclude that z = w which 
again contradicts our assumption. Therefore no such z exists and 
w is an atom of (M(A) , <;) . 

We next show that every atom of (M(A) , <;) is a column of W. 
Let a be an atom of (M(A) , ^) . Then a ^ and 

z^a => z=0 orz=a 

- 4. 
for z £ M(A) . By Lemma 5-8, a may be expressed as 



k 
a 



= > c. w 



(5-27) 



L l 
i=l 
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where c. is a non-negative integer, for i = 1, 2, ..., k. Since 
a ^ 0, there is some c. > 0, 1 ^ j ^ k, in this expression. 
It follows that w <, a. According to (5-27), w = or w = a. 
We know that w ^ 0. Thus a = w . 



Matrix W has at most A + 1 columns. This follows from the 
strict ordering, 



A = u^ > u 2 > • • • > u k = 



= v^ < v 2 < • • • < v = A 
Theorem 5-3 tells us that (M(A) , <) has at most A + 1 atoms. 



5.6 Generation of Normal Vectors and Vertices 



* 

We have seen that a normal vector y of I(P) at a vertex y 

is an atom of the partially ordered set (M(A) , ^) , where A is 

the cofactor matrix of A and Ax ^ b defines P. We have developed 

a method for generating the matrix W whose columns are the atoms 

of (M(A) , ^). In this section we complete our development of an 

integralization method for (2,2)-corner polyhedra by showing how 
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* 
all normal vectors and vertices of I(P) , hence of I(P), are generated, 

Our final product will be a system of inequalities A'x ^ b' which 

defines P* = H(I(P)). 

We begin by considering the following question. Given a 

i * 
boundary point y of I(P) , how do we tell which columns of W 

* i 
are normal vectors of I(P) at y ? 

i .* i v* 

Let y be a boundary point of I(P) . Since y € 1(P) we have 

■l 4 ' * ' 

y ■£. b. Let z = b-y . Then z ^ 0. z measures the "slack" 

c i 
of y . 

Recall that a normal vector y at y satisfies the following 

conditions : 



it) Y-(y-y X ) <: for all y £ I(P)' 
(ii) y 1 + Tv <= I(P)*. 



Using the slack vector z , we label two columns of W as shown 



below: 



A=u, >•••> u >"«>u >•••> u. =0 
1 p /a k 




= v 1 < ••• < v </..<• v < ... < v k = A 



v is the largest element in the lower row of W which is less than or 
equal to z.. u is the largest element in the upper row of W which is 
less than or equal to z„. 
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Theorem 5-4 



i * p q 

If y is a boundary point of I(P) then the columns w and w M 

i-1 i+1 * 

of W are normal vectors, y and y respectively, of I(P) at y . 



Proof 

* * 

According to Theorems 5-1 and 5-3, normal vectors of I(P) at y 

are columns of W. 

i+1 i 
We first observe that there exists a normal vector y at y 

i i+1 * 
such that y - T-y € I(P) • We shall consider all columns w of W 

i * i+1 

such that y - Tw £ I(P) to be candidates for y and show that 

a i \* 

all but w M fail to satisfy the condition w (y-y ) ^ for all y 6 I(P) . 

It will follow by elimination that y = w . 

For all columns w of W, y + Tw £ M(A). This follows from the 

following argument. Since w f_ M(A) , w = Ax for some integer vector x. 

Then Tw = TAx = ATx, so that Tw £ M(A). Since y 1 f M(A) , y X + Tw f M(A). 

Thus y 1 + Tw f. I(P) iff y 1 + Tw <; b. 

i i 1 

Now, for all columns w , where j < q, y -Tw J £, b because 



y„ + u. > b • For all columns \r , where q ^ j <; k, y - Tw ^ b 

because y. - v. <, b- and y_ + u. ^ b_. Thus y - Tw ^ I(P) 

for j = q, q+1, ..., k. However, for all column w , j = q+1, q+2, .-., k, 



w j • ((y 1 - Tw q ) -y 1 ) = w j .(-Tw q ) = 



V V. 

q j 
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because u > u. ^ and v. > v i 0. Thus all of these columns 
q J J q 

i i * 

fail to satisfy the condition vr > (y-y ) ^ for all y £ I(P) . 

Therefore, by elimination, y = */*. 

We next observe that there exists a normal vector y at y 

i i— 1 * 

such that y + 1y f I(P) . A proof similar to the one above 

establishes that y = w . HI 

In Theorem 5-4, the columns w^ and w need not be distinct. 
If they are distinct, then y is a vertex of I(P) . Otherwise y 1 

is merely a boundary point of I(P) . 

i * 

If y is a vertex of I(P) then Theorem 5-4 gives us two 

inequalities 

i-1 , i N _ 
ct . (x-x ) ^ 



i+1 ( i-\ n 

CT • (X-X ) £ 



(5-28) 



■{4-1 - »1 *4-1 

in the system A*x ^ b' defining H(I(P)). In (5-28), CT - = (A) y - 
and x = A y . 

We now consider the following question: given a vertex y of I(P) , 

how do we generate neighboring vertices, if any, of y ? 

i * ii 

Let y be a vertex of I(P) . We assume that z = b-y has been 

used to label two columns of W as before. We show how to generate 

slack vectors z = b-y and z = b-y of neighboring vertices 

i-2 , i+2 . , 
y and y respectively. 
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We refer to the previous labeling of W. If v > (that is , 



p > 1) then let 



i-2 i 

2 = Z - 



Tw P . 



(5-29) 



i-2 
Otherwise, disregard z If u > (that is, q < k) then let 



i+2 




i 




V 


z 




z 


+ 


u 

q 



TV 



(5-30) 



i+2 i-2 i+2 
Otherwise, disregard z . Using z and z , we label two 

additional columns of W as shown below: 




= v 1 < ... < 



... >u k 



= 



<v k =A 



v is the largest element in the lower row of W which is less than 
or equal to z^ . u is the largest element in the upper row of W 



which is less than or equal to z„ 



i+2 



Lemma 5-9 
i-2 



j i+2 *.- r t/t>%* 

y and y are vertices of I(P) . 
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Proof 

We prove that y is a vertex of I(P) by showing that w 

s * i+2 

and w are two distinct normal vectors of I(P) at y .We have 



i+2 i 

y = y 



V 



u 

q 



Tw<. 



Clearly y 1 <; b, so that y 1 f I(P) . Since w q «Tw q = 0, it follows 

q i+2 q i * 

that w • (y-y ) = w^» (y-y ) ^ for all y <= I(P) . Also, since 



l 
u 

q 



> 1, y 1+ + Tw q (i I(P) . Thus w q is a normal vector of 



* i+2 s * 

I(P) at y .By Theorem 5-4, w is a normal vector of I(P) at 



i+2 _. i+2 , s , q _ i+2 

q 



y . Since z„ < u^ we have w r w . Therefore y is a 



vertex of I(P) . 

A similar argument reveals that y is a vertex of I(P) . 



□ 



What we have described may be viewed as a second level 
two-dimensional division process in which the slack vectors 
of known vertices are the dividends, atoms of (M(A) , ^) are 
divisors, and slack vectors of new vertices are the remainders. 

In order to initiate the process we must have an initial 
vertex. This is easily obtained. In Section 5.1 we saw that 
I(tt,) D I(P) and I(tt 2 ) D I(P) are sets of boundary points of 
I(P), where tt, and tt 9 are boundary lines of half-planes defining P. 



133 



Furthermore we saw that I(tt, ) f] I(P) is the set of integer points 

11 1 
expressible as x -k..Ta , where x is an integer point in I(rr-) 

computed using the Euclidean algorithm and k 1 is an integer such 

that, 

b„-a •x 1 



b„-a -xt I 



We obtained a similar characterization of I(rr 9 ) ("| !(?)• 



If we take x -k-Ta , where k- = 



b„-a .x . .. 

, as being x , 



-J 



then we claim that x is a vertex of I(P). Consider y = Ax . 

11 11 

If we compute z = b-y we find that z- = and z„ < A« 

It follows from Theorem 5-4 that w and w are normal vectors 

s* 1 1 

of I(P) at y , where p = 1 and q > 1. Thus y is a vertex 

of I(P)*. 

Similarly, we can find an initial vertex drawn from 

I(tt 2 ) n i(p). 

1 1 * 

Using the slack vector z of the initial vertex y of I(P) , 

we label two columns of the matrix W as we did for Theorem 5-4. 

The circled slack vector z then points to the normal vectors 

2 ^1 
y and -y of I(P) at y . The division indicated by equation 

(5-30) is then performed on the labeled matrix W to obtain slack 

3 
vector z , which in turn is used to label a third column of W. 

3 
The circled slack vector z then points to the normal vectors 

2 4 
y and y — and so on. The diagram below illustrates the labeled 

matrix W. 
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... ( Y *-*) 



Example 5-6 Suppose W is given by 



W = 



73 56 39 22 5 3 1 
12 3 4 17 30 73 



and suppose z = f ). Then the labeled matrix W is shown below. 

X 70 ; 




Sample divisions are : 



(A) - (£H»J (-») 

(I) -(i)^J(-3)" 



Notice that Equation (5-29) enables us to check our division by 
working back from right to left. 



□ 
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The process terminates when a vertex y is found such that 
z~ =0, which must occur after at most k-1 divisions have been 
performed, where k is the number of columns of W. If the division 

process is executed from right to left, beginning with a slack 

-1 . -2n+l . 

vector z , then the process terminates when a vertex y is 

found such that z. =0. 

In the final labeled matrix W, each circled slack vector 

z points to the two normal vectors y and y of I(P) 

at vertex y . Each normal vector y is pointed to by the one 

* 
or two circled slack vectors representing vertices of I(P) 

at which y is a normal vector. 

The final result of all this is an alternating sequence 

2 4 2n 
Y y Y ••• Y 

1 3 2n-l 

y y ••• y 

•k 

of normal vectors and vertices of I(P) , which may be transformed 
into the alternating sequence 

2 4 2n 

O" 0" o" • • • X 

1 3 2n-l 

X X • • • X 

of normal vectors and vertices of I(P), where <y "~ =(A) Y — and 
x 1 = Ay 1 , for i = 1, 3, ..., 2n-l. We then write the following 
system of inequalities, 
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t 1 ^ n 

CT • (X-X ) £ 

ct • (x-x ) ^ 

ct 4 .(x-x 3 ) < (5-31) 



2n , 2n-l. 
a • (x-x ) «£ 0. 



1 9 1 

Since z.. = and z„ = 0, it follows from Theorem 5-4 that 

/AN , 2n /ON m ,- v -l 1 
Y = (^J and Y = (^J. Thus CT = (A) y = a and 

110 11 
ct .x = — Y *y = y- = b. , so that the first inequality 

in (5-31) is a^x ^ b^ Similarly a ln = (A)" 1 y 2n = a 2 and 

2n 2n-l 1 2n 2n-l 2n-l , , , 

ct .x = — y # y = y 2 = b„, so that the last 

2 
inequality in (5-31) is a «x ^ b„. These two inequalities 

are the two which define P. 

If we let Q denote the polyhedron defined by the inequalties 
(5-31), then it is easy to see that Q = H(I(P)). Since each of 
the half-planes whose intersection is Q contains I(P), it follows 
that I(P) s I(Q). Since the two half-planes whose intersection is 
P are included, it follows that I(Q) £ I(P). Thus I(Q) = I(P). 
Furthermore, by construction, every vertex of Q is an integer point. 
Thus Q is integral. Prom Lemma 1-1, Q = H(I(P)). 

The following lemma gives an upper bound on the number of 
vertices of H(I(P)). 
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Lemma 5-10 



H(I(P)) has at most A vertices, 



Proof 

We have previously argued that matrix W has at most A + 1 
columns. Therefore, in labeling the columns of W, at most A 
vertices are generated. 



□ 



We conclude this chapter with two examples illustrating 
the integralization of (2,2)-corner polyhedra. 

Example 5-7 Let P be defined by 



72-1 



[J a cy < [?] 



Then 



i - [.I f] 



["31 18" 
10 1_ 



w = 



T31 18 5 2 



Lo 



2 7 19 31_ 
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An initial vertex x of I(P) is given by, 



*° - » (1) - (A) 



where 



k = 



>-(D<?) 



31 



= -31 



Thus x° = (JJV and z° = ( 2 ° ? ). W is labeled as follows: 




The resulting system of inequalities defining H(I(P)) is 

© • KiS)) * ° 
(!) • HA)) * ° or 
ft) ■ (-(A)) * ° 
(1 ) • (*-©) - • 



2 


5 




X l 


£ 


72 


1 


3 




x 2 




41 





1 








11 


-1 


1 




4 


-5 


3 








2 



□ 
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Example 5-8 Let P be defined by 



488 216 
46 103 



279 
-317 



Using the Euclidean algorithm we find that gcd (488, 216) = 8 and 
gcd (46, 103) = 1. Since 8 { 279 the first inequality does not 
define an integral half-plane. The second one does. We replace 
the first by 61 x.. + 27 x„ £ I —3— I and obtain a new polyhedron 
P 1 defined by 



61 


27 




K l 




34 


46 


103 




x ? 


■£. 


-317 












L _ 1 



Then 



A = 



103 -46 
-27 61 



Matrix W has been computed in Example 5-5. An initial vertex x 







of I (P.) is given by, 



x 



- * (X) - (-V0 



where 



k = 



■»» -(x 4 0> 



X\) 



5041 



= 4 
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nn. 

Thus x = 



_ 2 6 8 2 ) and z° 



,^ on ). W is labeled as follows 
4/ol/ 




The resulting system of inequalities defining H(I(P..) is: 



61 
,27 



x- 



V38 



x- 



28 
■62 

4 
■8 



) • (*-Q) * • 



or 



61 


27 


9 


4 


2 


1 


17 


38 


46 


103 



34 

4 



-118 

-317 



□ 
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CHAPTER 6 



INTEGRALIZATION OF (2.2)-POLYHEDRA 



Let P be a (2, 2) -polyhedron defined by Ax <s b , where A is an 
integer m x 2 matrix of rank 2, and b is an integer m-vector. 
In this chapter we apply the integralization method of Chapter 5 
to the problem of integralizing P. 

We begin by reducing P to an intersection of integral half-planes 

i 1 i b ' 

by replacing each inequality a *x «s b. in Ax <; b by J a *x ^ 

i 
where d. = gcd (a.., a..). If we let t. denote the half-plane 

defined by a X »x «j b. and T-J denote the half-plane defined by 
b, 



d. 

1 



1_ i 

1. 

1 



d. a x * 



d. 

x 



then I Or.') = I(t.) for i = !> 2 > •••> m « Since 



mm m 

I(P) = I(T.) = fl I(t.'), P may be replaced by fl T.' without 
i=l X i=l L i=l L 

changing I(P). 

In the following we assume that this first integralization step 

has already been performed on P. That is, we assume that each row 

[a... a.„] of A has the property gcd (a.^ a ) = 1. 

Our plan is to integralize P by generating all vertices and 
normal vectors of I(P). Boundary points, vertices and normal vectors 
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of I(P) are defined just as they were for (2,2)-corner polyhedra 
in the previous chapter. We accomplish this by generating 
boundary points and normal vectors for a sequence of (2,2)-corner 
polyhedra which we shall call supporting corners of I(P). 

A (2, 2) -corner polyhedron Q is a supporting corner of I(P) 
iff I(P) c I(Q) and I(Q) and I(P) have one or more vertices 
in common. 

In the following we show how an initial supporting corner 
of I(P) is found. We then show how supporting corners of I(P) 
are used to generate vertices, normal vectors and other supporting 
corners of I(P). We conclude by showing that the process terminates 
after a finite number of steps, resulting in a system of inequalities 
A'x <z b 1 defining P'. 



6.1 Finding an Initial Supporting Corner 

Well-known methods associated with linear programming (see, for 
example, Dantzig [5]) may be used to find a vertex v of P. Associated 
with vertex v is a two-inequality subsystem Cx < d of the system 
Ax ^ b defining P, where v = C d. Cx ■$, d defines a (2,2)-corner 
polyhedron Q such that I(P) c I(Q). 





0" 


2 

a 




4 

a 


1 

X 




3 

X 
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Applying the methods of Chapter 5 to Q, we obtain the sequence 

2n 

a 

2n-l 

x 

of normal vectors and vertices of I(Q). Suppose x € I(P) for some i, 
where i = 1, 3, ..., 2n-l. Since a 1 - • (x-x 1 ) <; for all x f I(Q), 
it follows from the fact that I(P) c I(Q) that a 1- .(x-x 1 ) ^ for 
all x f I(P). Thus x 1 is a common vertex of both I(Q) and I(P), and Q 
is a supporting corner of I(P). 

Now suppose x 1 <£ I(P) for i = 1, 3, ..., 2n-l. In this case, 
P either does or does not intersect the boundary of Q 1 = H(I(Q)). 

Suppose P does not intersect the boundary of Q'. Then since 
the vertex v is in P and P is convex, it follows that P lies 
wholly within Q-Q' and thus I(P) = co and P' = cp , where cd denotes 
the empty set. 

Suppose on the other hand that P does intersect the boundary 
of Q'. Let e , for i = 1, 3, ..., 2n-3 , denote the edge 

2 i i -x / i+2 i N „ -s n -i 
[xPR x = x +X(x -x ) , «; \ <; 1 j 

of Q'. We argue that the intersection of P with the boundary 
of Q' is containing in exactly one edge e of Q', as illustrated 
in Fig. 6-1. 
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Fig. 6-1 



First, we see that neither the half-line edge 



{x € R 2 | x = x 1 + Mx^v), \ * 0} 



nor the half-line edge 



(x a 2 | x = x 211 " 1 + Mx 211 " 1 -v), \ * 0} 



of Q' contains any points in P. Otherwise, since v g P, the convexity 
of P implies that x £ I(P) or x 6 I(P) , contrary to assumption. 
Next we see that if two edges e. and e. of Q', where i < j , both 

contain points in P, then since v £ P, the convexity of P implies that 

i+2 i+4 j ^ ,./„% . - , 

x , x , ..., x € I(P), again contrary to assumption. Thus the 

intersection of P with the boundary of Q' is contained in exactly 

one edge e of Q 1 . We form a new (2,2)-polyhedron P by adding 

i+1 i 
the inequality a * (x-x ) «g to the system Ax ^ b and deleting 
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the subsystem Cx <; d. Since P-P lies wholly within Q-Q' and thus 
contains no integer points, we see that I(P ) = I(P). Thus we 
have reduced the number of inequalities defining P by one, 
without changing I(P). 

Computationally, these observations are implemented as follows. 

For each vertex x 1 of I(Q), i = 1, 3, ..., 2n-l, we compute the 

i i 

slack vector z X = b-Ax . If z ^ for some i, then x € I(P) and, 

as we have seen, x 1 is a common vertex of both I(Q) and I(P). Thus 

Q is a supporting corner of I(P). 

Otherwise z 1 ^ and x £ I(P) for i = 1, 3 , . . . , 2n-l. 

In this case, for each i, i = 1, 3, ..., 2n-3, we form the following 

system of in + 2 inequalities in X, given below in vector form as, 

z + Mz -z ) s: 

£ \ £ 1. 

(6-1) has a solution \ iff P intersects the boundary of Q' at edge 
e 1 . As we have seen, two cases are possible. In the first, (6-1) 
has no solution \, for i = 1, 3, . . . , 2n-3. In this case I(P) = 9. 
In the second case, (6-1) has a solution \ for exactly one i among 
1, 3, ..., 2n-3. We then form a new (2 ,2)-polyhedron P by adding 
the inequality o- 1+1 • (x-x 1 ) <; to Ax <; b and deleting the subsystem 
Cx «s d. 

If we have not found a supporting corner of I(P) or learned 
that I(P) =C P, we may apply the entire procedure described above 



(6-1) 
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to P . In general, during the k application of this procedure 
we either: (a) discover a supporting corner of I(P), or (b) find 
that I(P) = cp, or (c) obtain a (2 ,2)-polyhedron P defined by 
m-k inequalities such that I(P ) = I(P). 

We terminate this procedure when either (a) or (b) results 
from a given pass. If neither (a) nor (b) have resulted from 
the first m-2 passes, then we are left with a (2,2)-corner 
polyhedron P™ such that I(P m ~ ) = I(P). By definition, P™" 2 
is a supporting corner of I(P). 

We assume now and for the remainder of this chapter that 

th 
I(P) fi cp. Thus a supporting corner of I(P) is found on the k 



pass, where 1 ^ k ^ m-2. In the following, for convenience, 
refer to the resulting (2 ,2)-polyhedron P as P and to its 
defining system of inequalities as Ax ^ b. 



we 



6.2 Generation of Normal Vectors and Vertices 



k 
Let Q be a supporting corner of I(P) and let x be a common 

k k i 

vertex of both I(Q ) and I(P). We refer to the pair (Q , x ) as a 

state of the integralization process. Q is defined by two 

k k 
inequalities C x < d . We assume they are ordered such that 

det C > 0. We shall refer to the inequality toward whose boundary 



147 



line we are moving as the objective inequality . The objective 

inequality associated with Q is drawn from the system Ax <; b 

defining P. We shall denote the objective inequality associated 

. , „k , k , , . th . ., . . , 
with Q by a -x < b. , the j inequality m Ax < b. 

J k 
We proceed from state to state by means of a basic iterative 

step. A state which results from one application of this iterative 

step to (Q , x ), where the second (first) inequality of C x <, d 

is the objective inequality, is A + (-) successor state of (Q , x ) . 

The + successor of (Q , x ) takes one of the three forms illustrated 

in the diagram below. 




(Q\ x 1 ) 



k+1 i+2. 



(Q k± \ x 1 ) 



The + successor of (Q , x ) is characterized by a new vertex, 
a new supporting corner, or both. 



Example 6-1 Fig. 6-2 (a) and (b) illustrate a succession of states: 
(Q , x ) _> (Q , x ) 



(Q k+1 , x 1+2 ) 



.k+1 i+4 
(Q , x ). 



□ 



A state (Q , x ) used to start the process is an initial state . 
Q is defined by two inequalities C x ^ d both drawn from Ax <, b . 
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Thus either inequality may be designated as the objective inequality 

a «x <; b. . We saw how to find an initial state in the previous 
section. 



OBJECTIVE 
INeAUALlTH' 




OBJECTIVE ^ 
INEQUALITY 




(a) 



(b) 



Fig. 6-2 



A state which has no + successor or no - successor is a terminal 
state . 

We now describe the basic iterative step for determining the 

t k i 
+ successor of (Q , x ) , if one exists. 



^ For + successor, all instances of +(+) should be read as +(-). The 
opposite interpretation holds for - successor. 
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(1) Using the (2, 2) -corner polyhedron integralization method 
of Chapter 5, we compute the normal vector a ~ of I(Q ) at x . 

(2) We let 6X 1 = +T0- 1 - . We compute 

6y :L = A 6x X 
and compare §y with the slack vector z , where 

z = b-Ax = b-y . 

(3) If all components of 6y are non-positive then x + n 6x € I(P) 

i+1 i 

for all integers n ^ 0. a" is a normal vector of I(P) at x . We 

adjoin the inequality cr 1 - . (x-x 1 ) ^ to an evolving system of inequalities 

k i 
defining P'. The state (Q , x ) is a terminal state and has no + 

successor. 

(4) If one or more components of 8y are positive then we compute 
the non-negative rational numbers 



M = 






for all subscripts j such that 6y. > 0, and among these we identify 
a subscript $,,■< such that 



Xj k+1 * Xj 
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i t 
for all j such that 6y. > 0. We distinquish three possible cases 

corresponding to the three possible forms of the + successor of 

(Q\ x 1 ): (a) j k+1 = j k , (b) j k+1 ,1 j k and M k+1 * 1, (c) J k+] _ * j fc 



and <; M k+1 < 1. 



i+1 



(5a) (Ju + i = iO In this case \j, ^ 1. Since x + 6x £ I(P), 
is a normal vector of I(P) at x . We adjoin the inequality 



i+1 i 
0" — • (x-x ) ^ to an evolving system of inequalities defining P'. 



We let 



i+2 
x — = x 



1 + |>J 6X1< 



i+2 k k i+2 

x — is a vertex of both I(Q ) and I(P). We take (Q , x — ) as the 

+ successor of (Q , x ) . 



(5b) (J k+1 ^ J k and \J k+1 Ss 1) In this case, since 

x + 5x £ I(P), <j ~ is a normal vector of I(P) at x . We adjoin 

i+1 i 
the inequality <j ~ • (x-x ) <, to an evolving system of inequalities 

defining P'. We let 

i±2 i , ... . . i 

X = X + A.J 6x . 



[>j 



+ 

The following somewhat arbitrary rule insures the determinism of 

this step: If there is more than one subscript j such that \j 

is smallest, then choose j, if it is included — otherwise choose 

the smallest subscript. (Recall that j, identifies the objective 

k 
inequality associated with Q .) 



151 

k+1 
In addition, we let Q - be the (2,2)-corner polyhedron defined by 

a ~ • (x-x ) ^ and the new objective inequality a — «x ^ bj 

drawn from Ax <; b. We see that I(P) c I(Q ~ ) and that x - 

is a vertex of both I(Q - ) and I(P). Thus Q - is a supporting 

corner of I(P). We take (Q — , x — ) as the + successor of 

• ~k i* 
(Q , x ). 



( 5c ) k+1 * J k and ° ^ X Jk+l < ^ In thiS CaS6 ' since 

i i ~ k ±! 

x + 6x £ I(P), no new vertex is obtained. We let Q be 

the (2,2)-comer polyhedron defined by c ~ -(x-x ) £ and 

k+1 
the new objective inequality a _ «x «g b j k+ i dr awn from Ax <, b. 

k+1 i 

We see that I(P) e I(Q ~ ) and that x is a vertex of both 

I(Q - ) and I(P). We take (Q — , x 1 ") as the + successor of 
(Q k , x 1 ). 



This completes our description of the iterative step. 

The integralization process proceeds as follows. Starting 

with an initial state (Q , x ) , a sequence S of + successors 

k 2n 
is generated until a state (Q , x ) is obtained such that: 

(a) n ^ 2 and x = x 

or 

k 2n 

(b) (Q , x ) is a terminal state. 

If (Q , x n ) is a terminal state then starting with the initial state 

(Q , x ), a sequence S of - successors is generated until a terminal 

. -k* -2n\ . _ . , 
state (Q , x ) is obtained. 
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In case (a), the sequence S of states identifies n distinct 
vertices x , x , ..., x of I(P) and returns to the initial 
vertex x . If n > 2 then the generated inequalities cr i+1 . (x-x 1 ) ^ 0, 
for i = 0, 2, ..., 2n-2, define a bounded (2,2)-polyhedron P having 
n vertices. P is contained in P so that I(P) s I(P). Also, since 
each of these inequalities is satisfied for all x £ I(P), we have 
I(P) c I(P). Furthermore, each vertex of P is an integer point, 

so that P is integral. By Lemma 1-1, P = P*. If n = 2 then P' 

2 
is the line segment having x and x as endpoints. Only two 

inequalities are generated. They define the line containing P 1 . 

In case (b), the sequence S of states identifies n+1 distinct 

2 2n 9 

vertices x , x , ..., x of I(P). Beyond vertex x n we discover 

an infinite number of boundary points of I(P). It follows that 

P is unbounded. The sequence s" of states then identifies n'+l 

distinct vertices x , x , ..., x n of I(P). Beyond vertex 

-2n' 
x we discover an infinite number of boundary points of I(P). 

The generated inequalities, a •(x-x 1 ) < for i = 0, 2, . . . , 2n 

and CT .(x-x ) s for i = 0, -2 -2n', define an unbounded 

(2,2)-polyhedron P having n + n* +1 vertices. An argument 

similar to the one given for case (a) establishes that P = P 1 . 

In the following we show that the integralization process 

terminates after a finite number of states have been generated. 

Consider the sequence S of + successors of an initial state 

S + + + + 

(Q , x ). We divide S into subsequences S Q , S- , S_ , ..., where 
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S k * ^ or ^ = ^' *' ^' ***» consists of the consecutive block 

of states having the same supporting corner Q as first component. 

Each supporting corner Q is defined by a system of two inequalities 

k k k k 

C x ^ d . According to Lemma 5-10, (Q )* has at most A* = | det (Tl 

vertices. It follows that the subsequence S, contains at most 

A, states. 

Now consider the subsequences S Q , S- , S 2 , ... . Associated 

k 
with each subsequence is a supporting corner Q . Associated with 

corner Q is an objective inequality drawn from the system of m 

inequalities Ax ^ b defining P. This inequality is identified by 

the subscript j,. Consider the sequence j~, j-, j_ There 

are two cases to consider: 

(i) The sequence contains a j, such that j, = j. for some i < k. 
More simply, the sequence contains a repeat. We assume that j, is 

the first repeat in the sequence. In this case, it can be shown 

k 

that x is a vertex of I(Q ). It follows that, for some state 

(Q , x ) in S, , x = x . If n = then we have the degenerate 

case in which I(P) consists of the single integer point x . Otherwise, 

2 
n S: 2 (the case n = 1 is not possible since x is necessarily distinct 

from x ) and we have a state, as in case (a) above, at which S terminates, 

+ + 

Thus S, is the final subsequence of S . Since each member of the 

sequence j Q , j_, j 2 , ..., j has one of m possible values and j, is the 

first repeat, it follows that k ^ m. 
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(ii) The sequence j , j-, j», ... contains a j, which 
identifies an objective inequality drawn from Ax ^ b whose boundary 
line contains a half-line edge of P. In this case, some state 
(Q , x ) in S, is a terminal state, as in case (b) above. 
Thus S, is the final subsequence of S . The sequence j«, j- , . .., j 
contains no repeats, otherwise case (i) would apply. It follows 

that k < m. 

+ + + 
In either case, the number of subsequences S , S- , ..., S 

in S is at most m + 1. Thus S is a finite sequence of states. 

A similar argument can be made for a sequence S of - successors. 

It follows that the integralization process terminates after a finite 

number of states have been generated. 

We conclude this chapter with an example illustrating the 

integralization process for a (2, 2) -polyhedron. 



Example 6-2 Let P be defined by 
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We reduce P to the following intersection of integral half-planes, 
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An initial state (Q , x ) is found to be, 



Q°: 
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We generate a sequence of + successors of (Q , x ). The steps 

involved are displayed by giving a tabulation of the resulting 

i i 
z and Sy vectors, followed by the sequence of supporting corners 

and labeled W matrices. In the tabulation, we circle the component 

of z for which "k . is smallest. Inequalities defining P' are 

enclosed in boxes as they are generated. 
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The sequence of + successors terminates with (Q , x ) since 



16 
x = x . 



Fig. 6-3 illustrates P and H(I(P)). 
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contained in the non-negative orthant). Extending this algebraic 
characterization to include all totally integral polyhedra would 
be an interesting problem for future study. 

We have answered the second of the two questions above for 
all (n,k) -polyhedra, where k «- 2. We began in Chapter 4 by showing 
that an arbitrary (n,k) -polyhedron can be transformed into a 
(k,k)-polyhedron by a transformation which preserves integralization. 
This enables us to focus on the problem of integralizing (k,k) -polyhedra 
without loss of generality. We then described the solution to the 
simple case in which k = 1. 

In Chapter 5 we developed our main integralization results for 
the case k = 2. We considered the problem of integralizing a (2,2)-corner 
polyhedron. Using a transformation which preserves the hull-forming 
operation, we showed that vectors normal to the surface of the sought 
after integral polyhedron are necessarily atoms of a certain partially 
ordered set. We introduced a two-dimensional generalization of 
the division theorem for integers, and showed how this is used 
to generate the finite set of all such atoms. We then described 
a second level division process involving these atoms which 
generates all vertices and normal vectors of the integral polyhdron, 
and thus its defining system of inequalities. 

In Chapter 6 we considered the integralization problem for 
(2 ,2) -polyhedra in general. We showed how a sequence of supporting 
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(2,2)-corner polyhedra is generated and how the integralization 
method of Chapter 5 is applied to each, thereby producing all 
vertices and normal vectors of the sought after integral 
polyhedron, and its defining system of inequalities. 

The integralization process developed in Chapters 5 and 6 is 
number theoretic in nature. It makes use of simple arithmetic 
operations without recourse to exhaustive enumerative methods 
or search procedures. 

It had been our goal to generalize this integralization 
method to make it apply to (k,k) -polyhedra, where k is arbitrary. 
However, we encountered some difficulty in attempting to generalize 
to k = 3 which we now describe. 

Let P be a (3,3)-corner polyhedron defined by Ax «j b, where 
A is an integer 3 x 3 matrix with non-zero determinant and b is 
an integer 3- vector. We may integralize each of the three 

(3,l)-polyhedra (half-spaces) forming P. This gives us a 

1 11 

(3,3)-corner polyhedron P defined by a system A x «g b of 

three inequalities for which 

gcd [Ag] = gcd [Ag ; bg] 

for all face 1-subsets S. We may then integralize each of the three 

(3, 2) -polyhedra defined by two inequalities drawn from A x <: b . 

2 2 2 

Thus gives us a (3, 3) -polyhedron P defined by the system A x ^ b 

of three or more inequalities for which 
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gcd [A^] = gcd [Ag ;b^] 



for all face 2-subsets S. Thus far we have used one and two-dimensional 

2 
integralization methods to obtain P . The remaining integralization 

step is three-dimensional in nature and this is where the difficulty 

arises. 

To illustrate the difficulty in its simplest form, suppose P 

is a (3,3)-corner polyhedron defined by Ax ^ b, where 

gcd [A g ] = gcd [A g i b g ] 

1 2 
for all face 1-subsets and 2-subsets S. Then P = P = P . We regard 

* 
A as a non-singular linear transformation, and let I(P) be the image 

* 
of I(P) under A. Normal vectors and vertices of I(P) are defined 

using obvious generalizations of definitions of their two-dimensional 

counterparts. We let A be the cof actor matrix of A. As before 

it can be shown that normal vectors of I(P) are necessarily members 

of M(A) . However, it is unfortunately no longer true that normal 

vectors of I(P) are necessarily atoms of the partially ordered set 

(M(A) , £>. The following example serves as a counter-example. 



Example 7-1 Let P be defined by 

-4 11 
3-5 4 
3-3 2 



x 2 £ 3 
x 3 [ 3 
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It is easy to verify that 



gcd [A ] = gcd [A q |b ] 



S • S" 



for all 1-subsets and 2-subsets S. The cofactor matrix A is given by 



A = 



2 6 6 

-5 -11 -9 

9 19 17 



Fig. 7-1 below illustrates H(I(P) ), determined by actual construction. 
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The vector y = (4,2,2) is a normal vector of I(P) at vertices 

(3,1,3), (3,3,1) and (2,3,3). It is easily verified that 

Y £ M(A) . However it is also easily verified that (4,0,0) and 

(0,2,2) are members of M( A ) . Thus the normal vector y is not 

an atom of (M( A) , ^). _ 



Since normal vectors are not necessarily atoms, our 
two-dimensional integralization method does not generalize 
directly. In order to generalize the method it would first 
be necessary to identify some finite subset of M(A ) which 

contains all possible normal vectors of I(P) . Then a method 

** + 
for generating this finite subset of M(A ) would be required, 

along with a process for using this finite subset to generate 

•k 

all normal vectors and vertices of I(P) . In the two-dimensional 
case, it was the finite set of atoms of (M(A ) , g) that satisfied 
our first prerequisite. In the three-dimensional case, the 
appropriate finite subset of M( A) may possibly be those members 
of M( A) which occupy the first two ranks of the partially ordered 
set (M( A) , ^). This is merely a conjecture which seems to be 
a reasonable inductive inference to make from the two-dimensional 
case, once the set of atoms of (M(A ) , ^) are found to be insufficient, 

Our two-dimensional integralization method gives us reason to 
believe that there exist number-theoretic methods for integralizing 
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(k,k)-polyhedra for arbitrary k. Our experience in seeking a 
generalization leads us to believe that inventing these methods 
is apt to be difficult. The problem remains a fascinating one 
for us, and one whose solution we believe would be of great value 
in the field of combinatorial optimization. 
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